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We consider the corrections to the Boltzmann theory of electrical transport arising from the Coulomb inter- 
action in disordered conductors. In this article the theory is formulated in terms of quasiclassical Green's 
functions. We demonstrate that the formalism is equivalent to the conventional diagrammatic technique by 
deriving the well-known Altshuler-Aronov corrections to the conductivity. Compared to the conventional 
approach, the quasiclassical theory has the advantage of being closer to the Boltzmann theory, and also 
allows description of interaction effects in the transport across interfaces, as well as non-equilibrium phe- 
nomena in the same theoretical framework. As an example, by applying the Zaitsev boundary conditions 
which were originally developed for superconductors, we obtain the P(i5)-fheory of the Coulomb blockade 
in tunnel junctions. Furthermore we summarize recent results obtained for the non-equilibrium transport in 
thin films, wires and fully coherent conductors. 
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1 Introduction 

More than a hundred years ago Drude put forth a theory of metallic conduction QEl. In his theory 
electrons are treated as a gas of particles which are assumed to move along classical trajectories until they 
collide with one another or with the ions. The collisions abruptly alter the velocity of the electrons, as 
illustrated in Fig.^ Since that time several details of Drude's theory have been improved. In Drude's time, 
for example, it seemed reasonable to assume that the electronic velocity distribution at equilibrium was 
given by the Maxwell-Boltzmann distribution. About a quarter of a century later, Sommerfeld replaced 
the Maxwell-Boltzmann distribution by the Fermi-Dirac distribution. A further important step was the 
description of collisions beyond the relaxation time approximation, which finally led to the Boltzmann 
equation for the dynamics of the distribution function. In many cases the transport theory based on the 
Boltzmann equation successfully describes the electrical conductivity of metals. On the other hand the 
Boltzmann theory still assumes that electrons move along classical trajectories, and quantum mechanical 
interference effects are neglected. The latter become important at a high concentration of defects or in very 
small systems. Therefore, in the limits of strong disorder or small system size, the Boltzmann equation 
fails and the transport becomes non-classical. 

Alternatively, following Landauer [3 1, one can adopt a quantum mechanical description of charge trans- 
port: Consider, for example, a constriction as shown in Fig.|2] An incoming electron wave is transmitted 
or reflected at the constriction with amplitudes t and r. The current is then proportional to the transmission 
probability, i.e. the squared modulus of the amplitude. For simplicity, we restrict ourselves to zero temper- 
ature and a one-dimensional conducting channel for now. By assuming two different chemical potentials 
in the left and right reservoir, hl — (j,r + eV, the current through the constriction is given by 

I = 2e r L deAf(e)v{e)\t\ 2 ; (1) 
Jur 

the factor two is due to the electron spin. The product of the density of states Af(e) times the velocity v(e) 
does not depend on energy and is equal to Afv = 1 /2irh. As a result the conductance of the constriction is 
determined as 

G = 2^\t\ 2 , (2) 
with h/e 2 « 25.8 kfl. 

The connection between the scattering amplitudes and the conductance was first proposed by Lan- 
dauer 1 3 1. Generalizations of the conductance formula Q, including for example many transport channels 
and more than two leads, have been discussed intensively in the literature j4j |6| HI HI • The Landauer 
approach to the conductance has been successfully applied to describe transport through structures where 
the quantum mechanical coherence of the electron wave functions is maintained over the full system. A 
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Fig. 1 The classical picture for charge transport in metals: electrons move through the metal along classical trajecto- 
ries. From time to time they scatter and transfer momentum to the lattice. 




Fig. 2 The quantum mechanical picture for electron transport: an incoming electron wave is transmitted or reflected. 

fascinating example is the transport through quantum point contacts, which has been studied using both, 
semiconductor devices 110111 II as well as single atoms 1121 . 

In the following, we will focus on the transport in metallic systems, where the classical approach to the 
conductivity is still a good starting point, but quantum effects give rise to various types of corrections to the 
Drude conductivity. The two corrections that will be considered in this review are weak localization, which 
is a pure one-particle phenomenon, and effects arising from the combination of disorder and electron- 
electron interactions. 

In order to give a simple idea of how these quantum corrections arise and to highlight the connection 
between the semiclassical and Landauer's approach, we present here a physical picture of the origin of 
weak localization. For reviews on this subject see I13II14|[T51|16II17II18I . The weak localization correction 
to the conductivity arises as a result of the quantum interference of electron waves in disordered systems. To 
travel from point A to point B electrons move along classical paths. In order to obtain the total probability 
for a transfer from point A to B, in classical physics, one has to sum the probabilities for a particle to 
move along all possible paths. In quantum mechanics one has to sum the amplitudes of each path and take 
the modulus squared at the end. The relevant paths for weak localization are those with self-intersections 
and with the velocities of the incoming and outgoing paths in opposite directions as shown in Fig. [3] 
An electron can travel around the path clockwise or anti-clockwise. The two paths can be assigned an 
amplitude and the probability is then proportional to 

|*i + ^2\ 2 = l^il 2 + |* 2 | 2 + 2Rc(*^ 2 ). (3) 

The first two terms on the right hand side correspond to the classical probabilities. The third term, giving 
the interference between and $2, only appears in quantum mechanics. For the special type of paths 
considered, the interference term is always positive, i.e. due to interference the probability of the process 
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Fig. 3 A self-intersecting path as it is relevant for the weak localization correction to the conductivity. The dashed 
wavy line represents an inelastic scattering event which is responsible for phase breaking; t v is the phase coherence 
time, and D is the diffusion constant. 

is enhanced. Neglecting the interference corresponds to a classical description of the electrons (Drude- 
Sommerfeld theory, Boltzmann equation). 

This article summarizes recent contributions to the theoretical description of transport in mesoscopic 
conductors. The central questions addressed are the following: 

• In which way is transport modified by electron-electron interactions? 

• What happens in the presence of a large voltage, when linear response theory fails? 

Both questions are relevant for the interpretation of experiments: Shortly after the discovery of weak local- 
ization, it was found that similar effects in the conductivity are caused by the electron-electron interaction 
1 191 1201 12 II 1221 . Furthermore electron-electron interactions are also relevant for weak localization itself, 
since they provide a mechanism for phase breaking. However, while it is clear that inelastic scattering con- 
tributes to dephasing, the exact way this happens is far less obvious, as is evident from the recent debates 
in the literature 1 23l 1241 l25l 1261 l27l 1281 1291 1301 . The importance of studying non-linear transport is also 
apparent: Several transport experiments on mesoscopic samples are carried out at low temperatures. Un- 
der these conditions inelastic scattering events which drive the electron system towards local equilibrium 
freeze out. Therefore non-equilibrium conditions can be realized even at rather low voltages. 

Here we will use the quasiclassical formalism which, as its historical predecessor, i.e. the Boltzmann 
equation, is very powerful in dealing with situations far from thermal equilibrium. This was, for example, 
exploited in studying the dynamics of superconductors in the 1960s and 1970s 1 3 1 1 and hybrid mesoscopic 
systems in the 1990s 1 32 1. In particular, we will show, that the quasiclassical method provides a convenient 
framework to describe both quantum interference (weak localization) and Coulomb interaction effects in 
normal-conducting systems out of equilibrium. 

Outline 

The article is arranged as follows: In section |2] we introduce the Green's function technique in the non- 
equilibrium (Keldysh) formulation, which is the main theoretical tool used throughout this article. In 
particular, after giving the definition of the quasiclassical Green's function, we derive the equation of 
motion as well as the boundary conditions to be used in the presence of interfaces. Following the literature 
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we will demonstrate how a Boltzmann-like theory and the Drude conductivity are recovered within this 
formalism. In section [5] transport beyond the Drude-Boltzmann theory will be considered. By extending 
the formalism of section |2] contributions to the current density due to "maximally crossed diagrams" and 
due to the Coulomb interaction will be calculated. 

Sections |4]and[5] will be devoted to applications, i.e., the general expression for the current density of 
section|5]will be evaluated explicitly for different experimental setups. To begin with, in section|4] we will 
consider the limit of weak electric field and derive the well-known expression for the Coulomb correction 
to the Drude conductivity. This will be followed by a discussion of the the non-linear conductivity in films. 
Furthermore we will investigate the question of whether the phase coherence time t v , which is a central 
quantity for weak localization, is relevant for the Coulomb interaction corrections to the conductivity as 
well. In section|5]we will analyze non-linear transport in wires and in the presence of interfaces. In order 
to illustrate the power of the method we first show how to obtain the Coulomb blockade theory. We will 
also discuss the different regimes that arise depending on how the length of the wire compares with the 
thermal diffusion, inelastic electron-electron and electron-phonon lengths. Chapter[6]summarizes the main 
results. 

A few technical details are provided in the appendices. In appendix|A]we outline the correspondence 
of the quasiclassical formalism with the field-theoretic approach based on the non-linear sigma model. In 
appendix|S] we provide the technical details necessary to include the spin degrees of freedom. 



2 The classical theory of transport 

The traditional transport theory in metals is based on the Boltzmann equation. The central object is the 
distribution function /(k, r, t), where f(k, r, t)dr dk/(27r) 3 is the number of electrons (per spin direction) 
at time t in the phase space volume dr dk. The charge and current densities are given by 

/• Air 

P(r,t) = (-2e) J _/(k,r,t), (4) 

r Air P,lr 

JM) = (-2e) / ——f(k,r,t). (5) 

J (27r) J m 

In thermal equilibrium the distribution function reduces to the Fermi function. The Boltzmann equation 
determines the dynamics of the distribution function, and reads 

9 t / + v- V r /+iF-V k / = /[/]. (6) 

The left hand side of the above equation contains information on the energy spectrum, Hv = Vke(k), and 
the external forces F, whereas the collision term on the right hand side describes the scattering processes. 
For example, for impurity scattering, the collision term is given by 

'[/]=- / ^3^k,k'[/(k)-/(k')]. (7) 

In this section we will recall how a Boltzmann-like kinetic equation is derived within the Green's function 
formalism. To this end we use the non-equilibrium Green's function technique, as originally formulated 
by Keldysh 1331 . Our notation will mainly follow 1341 . After giving a few general definitions, we will 
introduce the quasiclassical Green's function. In the presence of impurity scattering, within the Born 
approximation for the self-energy, a Boltzmann-like kinetic equation and the Drude formula for the con- 
ductivity are recovered. To go beyond the semiclassical treatment and obtain the quantum corrections to 
the conductivity, to be discussed in sectionsp3l-|5j further approximations are necessary. From now on we 
set H = ks = 1) except in the final results, where, for clarity, we re-introduce the physical constants. 
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h 

Fig. 4 A closed time path contour, passing through t\ and t?.. 

2.1 The Keldysh formalism 

An object possessing a simple perturbation expansion is the so-called contour-ordered Green's function 
defined by 

G{x u x 2 ) = -i(T c ^( Xl )^(x 2 )), (8) 

where and vf^ are the usual fermion operators in the Heisenberg picture, and x» = (x*, ti). The brackets 
(...) denote an average over a statistical ensemble. The times ti vary on a (complex) contour which starts 
at some time to, and passes once through t\ and t 2 as shown in Fig.0] The contour-ordering operator T c 
orders the Fermi operators along the contour. Here we consider a special contour, the "Keldysh contour". 
It consists of two parts: The first part starts at to = — oo and ends at t = +oo, whereas the second part 
starts at +oo and goes back to — oo. The contour-ordered Green's function on this contour can then be 
mapped on a two by two matrix in "Keldysh space", 

G(x u x 2 )^G = ( g| ^ \ (9) 

with 

Gn(x u x 2 ) = -i(T*( Xl )&{x 2 )) (10) 

G l2 {x u x 2 ) = +i{^{x 2 )^{ Xl )) (11) 

G 2 i( Xl ,x 2 ) = -i{9(xi)^(x 2 )) (12) 

G 2 2(x u x 2 ) = -i(fy{ Xl )^{x 2 )}; (13) 

here T and T are the ordinary time-ordering and anti time-ordering operators. It follows from the definition 
that the components of G are not independent, and it is convenient to transform the matrix G as 



2V 1 1 / V _1 / V _1 1 

In this new representation the Green's function is 



QR qK 

G = ( qz qA ) } (15) 



with 



G R (x u x 2 ) = -i0(ti-t 2 )(*(a:i)*+(a; 2 ) + *t(a;2)*(a:i)) (16) 

G A (x\,x 2 ) = +ie(t 2 -t 1 )(^(x 1 )^(x 2 ) + ^(x 2 )^(x 1 )) (17) 

G K (x\,x 2 ) = -^(an)*^)-*^)*^)) ( i 8) 

G z (xi,x 2 ) = 0. (19) 
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Fig. 5 Graphical representation of the Green's function; we mark G with a shaded box, which can be understood 
as the Keldysh component of the self-energy, see Eq. J24l . 



To appreciate the physical meaning of the various Green's functions, let us consider first non-interacting 
electrons. The retarded and advanced Green's functions can then be expressed in terms of the eigenfunc- 
tions tp\ and eigenenergies e\ of the single-particle Hamilton operator, so that 

GfW(x 1) x 2 ) = V^M^l. (20) 
^ e - e x ± lO 

We emphasize that these two components of the matrix Green's function depend only on the energy spec- 
trum of the system, whereas the Keldysh component of the Green's function carries the information about 
the statistical occupation of the states. In thermal equilibrium the Keldysh component may be expressed in 
terms of the retarded and advanced components as Gf = [1 — 2/(e)] (Gf — Gf), where /(e) is the Fermi 
function. Generally, the equation of motion for G K constitutes the quantum-kinetic equation. Approxima- 
tions to this equation lead to the Boltzmann equation, Boltzmann-like equations, and generalizations. 
In order to find such an equation we start with the right-hand and left-hand Dyson equations which read 



(Go 1 - t)G = 6(xi - x 2 ), G(G 1 - E) = 8{x x - x 2 ), 



(21) 



where the Green's functions Go, G, and the self-energy E are considered as matrices in space, time, and 
the Keldysh index. Go is diagonal in the Keldysh index, the space and time dependence being given by 

(e = |e|) 



G 1 (x 1 ,x 2 ) 



1 2 
id tl - — (-iV Xl + eA(xi)) + e<f>(xi) + fx 
2m 



5(xi - x 2 ). 



(22) 



Disorder and interactions are contained in the self-energy. In the Keldysh space the self-energy has the 
same triangular matrix structure as the Green's function, 



E = 



E A 



Relation d23l allows us to express the Keldysh component of the Green's function as 

G K = G R Ti K G A . 
For a graphical representation see Fig. [5] 



(23) 



(24) 



2.2 The quasiclassical Green's function 

In this subsection we show how to obtain a Boltzmann-like equation. We start by defining the center-of- 
mass and relative variables 



x=i(xi+x 2 ), r = xi - x 2 



(25) 
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and Fourier transform with respect to the relative coordinate r in order to obtain the Green's function in 
the mixed representation 



G(xi, x 2 ) -> G(p, x) = / dre- ip - r G(x + r/2, x - r/2). 



(26) 



In the mixed representation, a matrix product of two objects A(xi,X2) and £?(xi,X2), as it appears for 
example in the Dyson equation, becomes 11351 



(A5)(p,x)=e 



A(p,x)B(p,x). 



(27) 



When both A and B are slowly varying in x a gradient expansion is justified, i.e. the exponential in the 
equation above is expanded and only the first term or first few terms are kept. 

By subtracting the right-hand from the left-hand Dyson equation and performing a gradient expansion, 
we obtain the equation of motion for the Green's function as 



id t + i— p • d_ 



G(p,x) = [£(p,x),G(p,x)] 



Here the scalar potential <fr and vector potential A are included in the gauge invariant derivatives, 

B t G = [9 tl +5 i2 -ie0(t 1 )+ie0(t 2 )]G tl)t2 (p,x) 

= [fy - i&j>(t + n/2) + ie<f>(t - ry/2)] G(p, rj; x, t), 

and 

9 X G = [9 X + ieA(ti) - ieA(ta)] G flt2 (p, x). 
In the equation above, t and 77 are center-of-mass and the relative time, defined by 

!, 

We introduce next the ^-integrated (quasiclassical) Green's function 



5tlt2 (p,x) - iy dedre- 1 P r G(x+^,ti;x 



g(p,?7;x,t), 



(28) 



(29) 
(30) 



(31) 
(32) 

(33) 
(34) 



where ^ = p 2 /2m — fi, and p is a unit vector along the momentum. The ^-integration is understood as 
a principal value integration, which implies that occasionally suitable (equilibrium) quantities have to be 
subtracted, in order to ensure convergence. In the entire article we will keep the notation of small g for the 
^-integrated Green's functions, and capital G for the original Green's functions. When approximating the 
density of states as an energy independent constant, the ^-integration is related to the original integration 
over the momentum p according to 



dp 

(2tt) 3 



iV / d£ 



dp 

47T ' 



(35) 



Neglecting the ^-dependence of the self-energy E, and integrating (12 8t leads to the equation of motion for 
the quasiclassical Green's function (Eilenberger equation |36|) 



d t + v F p ■ <9 X 



,g(p,x) = -i [E(p,x),.g(p,x)] . 



(36) 
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We note that the non-homogeneous term appearing on the right-hand side of the Dyson equation has can- 
celed when deriving the equation of motion for the quasiclassical Green's function. As a consequence, the 
latter is determined only by a multiplicative constant. Therefore the equation of motion has to be supple- 
mented by a normalization condition for the Green's function, which turns out to be of the form "g 2 — 1", 



J dii<? ttl (p,x)g tlt ,(p,x) = 8{t-t')l, 



(37) 



where 1 denotes the unit matrix in the Keldysh space 1 36 37 1. We now give some relations that are specific 
for impurity scattering, when it is treated within the self-consistent Born approximation. By assuming a 
Gaussian, (^-correlated impurity potential one observes that the impurity self-energy is related to the s-wave 
part of the quasiclassical Green's function, 



S tl t 2 (P, x ) = -- 



dp 

2r / 4?r 



9tit 2 (p,x), 



(38) 



where r is the usual scattering time. Due to the matrix structure of the theory, the above equation is satisfied 
by all three components (retarded, advanced, Keldysh) of the self-energy. The retarded and advanced self- 
energies have a simple structure, S^j = =F(i/2r)5(ii — t 2 ), and the retarded and advanced Green's 
functions are given by 



9 5f)( P! x)=±^ 1 -t 2 ). 



(39) 



We already mentioned that the equation of motion for the Keldysh component of the Green's function 
constitutes the kinetic equation. By also transforming to the mixed representation for the frequency/time 
variables the kinetic equation reads in the present case 



d t +v F p-d x g K (p,e;x,t) 



<, h [p.<:X.I) - / ^V(p, e;x ,i) 



(40) 



which clearly reminds us of the Boltzmann equation. A direct connection with the Boltzmann equation has 
been suggested in 1 38 1. By defining the distribution function as 



/(p,x,i) = 



(41) 



and ignoring the explicit dependence of g on p so that the derivatives of / are approximately given by 



V P / = 



V x / = -\[V^g K -e{V^)d,g K ) 
dtf = -\[d t g K - e{d t m9% 



(42) 
(43) 
(44) 



one recovers from the kinetic equation J40l > the Boltzmann equation (jfji with the external force given by 
F = e(V x + d t A). 

Besides the kinetic equation we also need a rule to calculate the physical observables. Generally, the 
charge and current densities are related to a combination of all components of the Green's function, 



p(x,t) = ie [G K (x,i;x,i) - G fl (x,i;x,i) + G A (x, i; x, £)] 

j(x,t) = [V x — V X ' +2ieA(x, t)] [G K — G R + G A ] (x, t; x', 
2m 



(45) 
(46) 
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In terms of the quasiclassical Green's functions, the charge and current densities are 1 34 



p(x, t) = 2e7V, 



ir f de f dp K , A . . . 

2 J 2^ J l^ 9 (p^;*.*)-^*) 



(47) 



de f d P„. - JO 



j(x,i) = &ttNq j — J ^^ P3 K ( P! e;x,i). (48) 
2.3 Diffusive limit 

In this paper we will focus our discussion on the diffusive regime, characterized by small external frequen- 
cies and momenta such that ujt, qvpr <C 1. Under these conditions the Green's function becomes almost 
isotropic and may be expanded in spherical harmonics. In the present case, it is sufficient to retain the s- 
and p-wave components, 

ff(p,x) = &(x)+p-gp(x). (49) 

By inserting Eq. ( I49i into the kinetic equation and separating the s- and p-wave parts, one obtains the 
relation 

g p (x) = -v F rd x g s (x), (50) 
and finally the kinetic equation for the s-wave component of the Green's function 

(5 t -Dd x d x )g?(e;x ) t)=0 i (51) 

where we have introduced the diffusion constant D = v^t/3. From Eqs. ( I47i - ( 1501 . an explicit expression 
for the current density follows: 

j(x, t) = -enDN d x I ^-gf (e; x, t) = -i^V x p(x, t) + 2e 2 ^7V E(x, t), (52) 

in agreement with the Drude formula for the conductivity. 



2.4 Boundary conditions 

In particular, we will be interested in the transport properties of finite systems, i.e. metallic systems with 
boundaries to insulators and reservoirs. We note that, due to the ^-integration and the gradient expansion, 
the quasiclassical equations are not directly applicable across sharp interfaces and therefore must be sup- 
plemented by boundary conditions. For a surface with specular scattering, and for the non-interacting case 
we consider at present, the matching conditions read 

<?(p;x_) = i?<?(p;x_) +T5(p;x+), (53) 
g(f>;x+) = i?.g(p;x + )+7g(£;x_), (54) 

where T and R are the transmission and reflection probabilities at the surface, p and p are the directions 
for incoming, transmitted, and reflected particles as shown in Fig.|6] From the sum of the two equations 
and using the relation R + T = 1 it is found that the antisymmetric combination of g(p) and g(p) is 
continuous at the interface, 

<Kp; x -) -<Kp; x -) =5(p; x +)-<7(p;x+), (55) 

' Note that we use "T" for both, the transmission probability and the temperature. The meaning will be apparent from the 
context. 
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R, p 




Fig. 6 An interface with specular scattering. An incoming particle moving in direction p is reflected with probability 
R into p, and transmitted with probability T into p; 6 is the angle between the trajectory and the normal vector n of 
the interface. 



from which one observes that the current normal to the interface is conserved, since 



^;(n-p)5(p;x_ ; 



^(n-p)g(p;x+). 



(56) 



After subtracting J54i from (I53i and making use of the continuity of the antisymmetric combination of the 
Green's functions it is seen that the symmetric combination of the Green's functions jumps at the interface, 
with the size of the jump determined by 

<?(p;x-) -<?(P> X -) = 2^ { <?(P; x -) +5(P,x-) - 5(p;x+) +5(p,x + ) | . (57) 

The generalization of this boundary condition to the interacting case is given below, compare Eq. (1148b . 

After these general considerations we now focus our attention on the diffusive limit. Equation \51\ 
is valid both for clean and dirty metals. In a dirty metal, however, a matching condition which involves 
only the angular averaged Green's function g s (x) instead of g(p;x) would be desirable. Recall that in 
the diffusive limit the angular dependence of g(p, x) can be taken into account by keeping only the s- and 
p-wave parts. One might then be tempted to expand the Green' functions in (1571 in an s-wave and a p-wave 
part, multiply d57l > by p, take the angular average and hence obtain the matching conditions in the diffusive 
limit in terms of the s-wave andp-wave components of the Green's functions. However, the task of finding 
the matching condition turns out to be not that simple, since close to the interface the angular dependence 
of the transmission and reflection coefficients may generate higher harmonics in the Green's functions, and 
it is not justifiable to terminate the expansion after the first two terms. Only at some distance from the 
interface will such an expansion hold. The task is thus to derive an effective condition which connects the 
Green's functions "far" from the interface, i.e. on distances of the order of the elastic mean free path. By 
integrating the kinetic equation and making use of \51\ it has been shown 1391 1401 14T1 that this matching 
condition can be written as 



-D(h ■ <9 x )<?f (x) = r [g? (x_) - gf (x+)] , 



(58) 



where T is a function of the angular dependent transmission and reflection amplitudes R and T. For a 
strong barrier (T/ R <C 1), T is given by 



vp I u(cosfl) C0S# ; 

/ V ' 2i?(cos< 



while for a weak barrier (R/T <C 1) 



[ d(cos0) cos 3 (0)^ 
Jo T 



(59) 



(60) 



12 



P. Schwab and R. Raimondi: Charge transport in disordered interacting electron systems 



see 1 41 1; for a barrier of arbitrary strength the parameter T may be obtained from the solution of an integral 
equation ATI . 

As a simple application of the matching condition J58i we consider a system consisting of a diffusive 
wire of length L which is attached to two reservoirs. We study the system in a non-equilibrium situation 
with an applied voltage Vi — V r = V, where the subscripts / and r indicate the left and right reservoirs, 
respectively. The classical resistance of the structure is the sum of the wire resistance and the interface 
resistances, i? to t = -Rwh-c + Ri + R r , and the current is / = V r /i? to t- We will demonstrate now how this 
result can be reproduced within the quasiclassical formalism. 

Within the quasiclassical formalism the current flowing in the wire is given by 

7 wiro = -enDNoA ^d x g*(e,x), (61) 

where A the cross section, and x = ... L is the position along the wire. The reservoirs are assumed to be 
in thermal equilibrium, with the distribution function given by gf' l ' r (e) = 2 tanh [(e + |e|V; !r )/2T]. The 
boundary conditions read 

-Dd x gK(e,x)\ x=0 = iMsf'^-sfM)] (62) 

-Dd x g?(e,x)\ x=L = T r [g?(e,L)-g«> r (e)}, (63) 

where Ti, r are the transparencies of the left and right interface. The current across the interfaces can be 
written as 

Il = -i e r^7V |d e [ 5 f( e ,0)- 3 f ; (e)] (64) 

Ir = -\eT r AN Q Ue[gf' r {e)-g?{e,L)]. (65) 

Furthermore we assume that the wire is so short that we can neglect inelastic scattering. In this case the 
kinetic equation for the distribution function inside the wire simply becomes —Dd x gjf = 0. Solving this 
equation one finally finds the current as / = V/(Ri + R W i IC + R r ), with 

Rr = l/(2e 2 .4JVor0 (66) 

i? wirc = L/(2e 2 AN D) (67) 

Ri = l/(2e 2 ^iV r r ) (68) 

as one would expect for three resistors in series. 



3 Quantum interference in diffusive conductors - general formalism 

Quantum effects give rise to deviations from the classical expression of the electrical conductivity of a 
metal: The conductivity depends on the sample specific realization of the impurity potential, and even after 
averaging over all possible realizations of the impurity potential corrections to the conductivity remain. 
The quantum corrections to the average conductivity in a metal with diffusive electron motion are weak 
localization, and the interaction contributions to the conductivity. The latter are often classified as the 
particle-particle (Cooper) channel, which is related to the exchange of superconducting fluctuations, and 
the particle-hole channel, which is related to the exchange of charge fluctuations (spin singlet channel) or 
spin fluctuations (spin triplet channel). In this article we concentrate on the average conductivity, and in 
particular on the Coulomb interaction in the particle-hole channel. We will neglect the Cooper channel. 
This is justified in non-superconducting metals, since in this situation the relevant interaction parameter 
scales down under the renormalization group flow. For completeness we will also discuss briefly weak 
localization. In this section we will give the general expressions for the contribution to the current density 
due to weak localization and due to the electron-electron interaction. The following sections |4] and [5] will 
contain specific applications. 
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3.1 Weak localization 

In a weakly disordered metal, quantum interference leads at low temperature to deviations from the Drude- 
Boltzmann theory. Gorkov et al. | 42 1 and Abrahams et al. | 43 1 showed that the summation of maximally 
crossed diagrams gives rise to divergences in the conductivity for arbitrarily weak disorder in dimensions 
less than and equal to two. This so-called weak localization correction to the conductivity is due to electrons 
diffusing along closed paths, where quantum interference causes an enhanced backscattering probability, 
as discussed in the introduction. The weak localization contribution to the current density is given by 



4 f°° 

SiwL(t) = -e 2 D- J d v (x, x)E(t - 77), 



(69) 



where D is the diffusion constant, the short time cut-off r in the ^-integration is the elastic scattering time, 
and C* _„(x, x) is the cooperon at two coinciding points in space. In the presence of a vector potential, 
the cooperon is given by the solution of the differential equation 

2^ - D(V» + ieA c ) 2 j C* )V (x, x') 

= 5(x-x')5(r)-r)'), (70) 

with Ac = A(x, t + 77/2) + A(x, t — rj/2). We recall now the results for the conductivity. In the absence 
of external fields the cooperon at two coinciding points in space is given by 

where d is the spatial dimension and a phenomenological dephasing time t v has been introduced; micro- 
scopically dephasing arises from inelastic scattering. Inserting the above into Eq. ( I69K one arrives at the 
standard result O El El El 



-hV^P (d=l) 
<Wl = <( -^nHr v M (d = 2) (72) 



2ir' 2 h 



y/l/Dr\p + const (d = 3). 



The connection to the Green's function formalism is the following: In section|2]we demonstrated that, by 
approximating the impurity self-energy by Xsom, the Boltzmann equation for the distribution function and 
the Drude conductivity are recovered. The kinetic equation including the weak localization correction has 
been derived in |44|, considering also the maximally crossed diagrams, 



with 



Smc, (73) 



(74) 



E mc - /,-' /\ \ + /--'/<- ^ ^ + //V'.'^^- + • • • (75) 

The self-energy has to be determined self-consistently. The contributions of the crossed diagrams to the 
retarded and the advanced self-energies are small and are therefore neglected |45|. The contribution to 
from the diagram with m impurity lines contains the sum of products with (2m — 1) Green's func- 
tions, Xm=r7 (G R ) n G K (G A ) 2m ~ 2 ~ n , i.e. integrals involving all combinations of the Green's function 
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Fig. 7 Relevant graphs for the Keldysh component of the self-energy, E^ c ("mc": maximally crossed). The first 
graph leads to the term proportional to g K (— p) in the kinetic equation, the other two graphs are responsible for the 
term containing an angular average over g (p). 



from G K G A . . . G A to G R . . . G R G K . The momentum integrals reduce to integrals over pairs of Green's 
functions of the form 

and the analogously defined integrals r\ AA and r\ RA . For small q the integrals rj RR and r\ AA are of the 
order 1 / (ej?r), whereas the integrals 7] RA are of order one, 

T) RA fa 1 - t[-UJ + Dq 2 ]. (77) 

Terms involving i] RR or 7] AA will thus be neglected. Following this rule the graphs contributing to E mc are 
shown in Fig.0 In order to understand how the above diagrams modify the kinetic equation, let us examine 
the first diagram in Fig. Q In this diagram there is an equal number of retarded and advanced Green's 
functions and a Keldysh Green's function in the middle, i.e., we have a sequence G^ . . .G^ G^ G^ ...G^ . 
The self-energy of the first diagram of Fig.0has the form 

5E K (p, e; q, W ) - C(Q, uj)G k (Q - p, e; q, u), (78) 

where Fourier variables p, e and q, uj refer to relative and center-of-mass coordinates, respectively, and the 
integration over the cooperon takes care of all the series of maximally crossed diagrams: 

C(Q,u) = T [r, RA + (r, RA ) 2 + ...]=rr, RA /[l-r, RA ] (79) 

- ^TW (80) 

To finally derive the modified Eilenberger equation, the commutator between the self-energy and Green's 
function, which appears on the right-hand side of Eq. d28i . has to be integrated with respect to £. As a 
result we have 

\ J d^£,G]* = i^^£C(Q,u,) / dtG«(G R -G A ). (81) 

Q 

Following ( 44 1 we use the Ansatz 

G K (p,q) = iff A '(p,q) (G R (p + q/2) - G A (p- q/2)) , (82) 
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which is based both on the fact that G R and G A are not corrected beyond the Born approximation, and that 
the equation above gives the correct quasiclassical Green's function upon ^-integration. A similar analysis 
can be carried out for the other diagrams of Fig. after some algebra the kinetic equation is found as 



d t +v F p-B x 3 K (p,ry;x,t) 



9 (p,??;x,*) 



4tt 



9 OMjx,*) 



a(t,t') = 



dt'a(t, t' 
2 



5 K (-p,ry;x, t') 



4tt 



c. 



t-rf/2 



,(X,X), T)' = t-t', 



(83) 
(84) 



where we put the general version of the cooperon ^, (x, x) with three different time arguments as given 
in J70l >. The maximally crossed diagrams give rise to an additional scattering term on the right hand side 
of Eq. J83i which is non-local in time, with the dephasing time as the relevant time scale. 

We restrict ourselfes now to the diffusive limit and calculate the change in the current density due 
to the maximally crossed diagrams. One identifies two contributions. The first term appears because the 
additional scattering term may change the solution of the kinetic equation for the s-wave part of the Green's 
function, 



-,K 



8gf , and 



<5jBorn(x,£) = -enDAf d x / tt 5 9s (e; x ,*) 

J 27T 

= -.DV x <fy)(x,i). 



(85) 
(86) 



The second term arises from the modified relation between the s-wave and p-wave contribution to the 
Green's function, 



g p = -VFTd x g s + vft 
which leads to 

$jmc(x,t) 



dt'a(t,t')d x g s 



2eDr J d^C^lf (x, x)9 x J —g* (e; ^t-rf) 
— J d77C^ /2 (x, x) [- V x /j(x, t-rj) + 2e 2 7V E(x, t - r,)] 



ttN ( 

The total weak localization correction to the current density is thus 

^JWL = <5jBorn + <5jmc- 



(87) 

(88) 
(89) 

(90) 



In this equation the sum of both terms is needed in order to ensure charge conservation. In the special case 
where the electric field and the charge density are homogeneous in space, the weak localization correction 
to the current density, as given in Eq. ( I69l l. is recovered. 



3.2 Interaction correction to diffusive transport 

Shortly after the discovery of weak localization it was found I20l l21ll22l that similar effects in the conduc- 
tivity are also caused by the electron-electron interaction. The interaction correction to the conductivity in 
the particle-hole singlet channel, for example, is given by [27 1 



<5oeei 



4e 2 L> 
nd 



drj 



siuh^Tr]) J \AirDr] ) 



■nTrj 



d/2 



(91) 
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which leads to 




-4M^j-^/hD/k B T (d 



1) 
2) 
3) 



(92) 



1.83^^^ k B T/hD + const, (d 



The inclusion of the triplet channels does not change the functional form of the temperature dependence 
of the correction to the conductivity, but modifies the prefactor, which then depends also on the strength of 
the electron-electron interaction in the spin triplet channel |46 47 1, see also Eq. J122i below. 

Whereas a simple and convincing physical interpretation of weak localization exists, we are not aware 
of as simple an interpretation of the interaction effect. However, attempts have been made [48, 49 1, and 
we present the main ideas. First, one observes that the impurities perturb the charge distribution in the 
metal, n — > n + Sn. For example, a single impurity gives rise to the so-called Friedel oscillations of the 
density, which persist even at large distance Sn(r) ~ am(2kFr)/r 3 . In the presence of many impurities 
the electron density becomes non-uniform, with the details depending on the particular distribution of the 
impurity positions. When both electron interaction and disorder are present, the charge inhomogeneity acts 
as an additional scattering potential, which, within the Hartree approximation, is given by 



It is clear that this additional scattering potential may affect the elastic mean free path. The way this hap- 
pens, however, is less obvious. In the following we give an argument 1 4-8 1 which shows that the interaction 
contribution to the conductivity is due to quantum interference. Let us consider two classical paths an elec- 
tron can follow to travel from point A to point B. Let path "1" and "2" be identical up to an extra closed 
loop in path "2", so that there is a phase difference between the two amplitudes, "J 2 = exp(i(/?i oop )\['i . The 
sign of the interference term Re^*^) is positive or negative, depending on the phase ^i 00 p, and there- 
fore processes of this type, when averaged over all the impurity configurations, give a negligible small 
contribution to the total probability for traveling from A to B. This will be different in the presence of 
electron-electron interactions. To first order in the Hartree potential Vh we may write the interference 
term as Rc^^^Vh-). Now note that the charge inhomogeneity 5n(r), and therefore Vh (r), is related to 
(virtual) electrons or holes which propagate along closed paths, since 



In the case where a virtual hole goes around the same closed loop as the path number "2", the phase factor 
exp(i</?i oop ) cancels, and ^1 and ^2 interfere coherently. A graphical representation of the relevant process 
is shown in Fig. [8] 

The Coulomb correction to the electrical conductivity was investigated in the past with various methods, 
ranging from standard many-body diagrammatic calculations of the Kubo formula to a field-theoretic de- 
scription based on a mapping of the original problem of interacting Fermions to a non-linear sigma model 
for matrix fields, and also to a generalized kinetic equation approach [ 19 46 47lll3ll34ll50lll8M51ll52ll53l . 
In the following we will briefly outline how the Coulomb interaction in a disordered medium can be incor- 
porated in the quasiclassical formalism we employ in this article. We start from the kinetic equation and 
the expression for the current density in the diffusive limit, 




(93) 




(94) 



d t g s - Dd x (g s d x g s ) = 



(95) 



and 



j(x,t) 



eirDN (g s d x g s ) K {x,t). 



(96) 
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Fig. 8 Graphical representation of processes leading to the interaction correction to the conductivity. The figure shows 
two paths for the propagation from A to B. One of the paths contains an extra closed loop. The corresponding phase 
factor exp(iyji 00 p) cancels with a phase factor due to a virtual particle which is created at (xi,ti — rj), propagates 
along the same closed loop, and is finally absorbed at (xi, ti). 

These two equations generalize the expressions we gave in section|2]to the situation where the quasiclassi- 
cal retarded and advanced Green's functions are not necessarily equal to plus or minus one. The Coulomb 
interaction is now introduced by adding internal electric fields to the external one, 

0^ </> oxt l + int , (97) 

and physical observables are obtained after averaging over the internal fields. Formally this can be carried 
out by means of a Hubbard-Stratonovich transformation within a functional integral formalism. More 
physically, this approach corresponds to describing the Coulomb interaction in terms of the exchange of 
scalar and longitudinal photons. The internal field has a non-trivial structure in Keldysh space, 



lint _ I 91 92 
t>2 9l 



(98) 



and the field fluctuations are related to the retarded, advanced and Keldysh components of the (screened) 
Coulomb interaction, 



"ie 2 ( )™( I = ^ ( \,a ' n ) • (99) 



l) (0102) \ _ 1 ( V K V R 



0291) \9292) 



2 V V 



In this language, screening effects appear as self-energy corrections of the photon propagators. To stress the 
difference between an external classical field and an internal quantum one, we observe that <f>i and cf>2 are 
related to the sum and difference of the values that the field takes in the upper and lower parts of the Keldysh 
contour, respectively. An external field takes the same value on both parts of the Keldysh contour and hence 
its (12)- and (21)-components in Keldysh space vanish. Finally we define the distribution function F via 
the s-wave part of the quasiclassical Green's function according to 

gf = gfF - Fgf. (100) 

We now show how to calculate quantities like the Green' functions, the density of states or the current 
density perturbatively in the interaction. For the Green's function, for example, we write 

sd*""]) = f J -°) + ($ £). 
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Note that due to </> 2 ^ 0, the (21)-component, i.e. 5g z is nonzero. Causality is restored only after averaging 
over the internal field fluctuations, and the (21)-componentthen vanishes. The normalization condition for 
the Green's function, g s g s = 1, allows us now to express the variation of the retarded and advanced Green's 
function and its average over the internal field as 



5g R 
Sg A 



Z K . 



(5g A ) q 



\{5g K 8g z )^, 
l -{5g z 5g K )^. 



(102) 
(103) 



In order to determine the correction to the retarded Green's function to first order in the screened Coulomb 
interaction it is sufficient to determine 8g K and Sg z to first order in the internal field. The product of Sg K 
and Sg z also enters the current density, since 



-enDN Q ((g s d x g s ) K }^ 



-enDNo 



d x (2F + (Sg K )^) 



+ (Sg K (d x Sg z ))^ t F a + F a ((d x 6g z )6g K )^ 
To first order in the field <p mt one finds from the equation of motion 

(d t + DdJ*)Sg z (ri;x,t) = 2i<fe(x, t)S(t]) 



(d t - DdJ y 



(104) 
(105) 

(106) 
(107) 



where the time derivative d t contains the external field <j> cxt , but not the internal field cj> mt . The ". . ." 
indicate a term proportional to 2 ; this term is of no importance here, since the correlator (0 2 (j>2) is zero - 
The times r\ and t are the relative and the center-of-mass times. In order to solve the above equations it is 
useful to define the diffusion propagator (the diffuson) via the following equation: 



[dt - DdJ x } D"(x, t; x', t') = <5(x - x')6(t - t'), 
which allows us to express the product of 5g K and Sg z as 



dt'(^(x)^ ti ( Xl )) 0i , 



(^tti( x )^ti( x i))*'"' 
2i J dx'dx"dt'dt"D t - tl 



(108) 



(109) 



xF I t - ti;x',t' 



t - U 



t + ti , , t — t] 
-2-;x,t - — 

v R (x',t'-y,t") 



(HO) 



As an immediate application of this equation one may obtain the correction to the density of states in 
equilibrium. The translational invariance both in time and space allows simplification of the above multiple 
integrals, and one arrives at 



SN(e) 



= Re((5g K (e)) 01 nt = Rc 



El 'tt- tanh 
J 2tt 



2T 



1 

-iu + Dq 2 



V R (q,uj) 



■ (HI) 
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By collecting the various pieces, we find for the current density in the presence of Coulomb interaction 
the result <5j' E ei = Sj Bo m + Sj v , with 

SjBor n (x,i) = -enDN V x J ^-(5g K (e;x,t))^ t (112) 

5jV(x,t) = e4irDN J dyy dx' At' dx" At" 

xRe[fb(-»?; x, t ~ r//2) J D r '(x, t - V /2; x', t' - r,/2)F Q { m x', t' - V /2) 
xV K (x', t>- x", t")(-iV x ) J D°(x", t"; x, t - 77)] - (113) 

Notice that this expression of the current density is valid for an arbitrary form of the distribution func- 
tion and of the diffuson. This allows us to investigate electrical transport in different experimental and 
geometrical setups as will be carried out in the following two sections. 

The expression for the current density J 1 1 31 has been derived first in [54|, both using diagrammatic 
techniques as well as with the Keldysh version l51M52l[T25l of the non-linear sigma model [46 1. Equation 
J 1 1 31 generalizes earlier results, which are valid in the absence of the external vector potential [56) or 
near local equilibrium f27ll55l . In |55| the spin-triplet channels and Fermi liquid renormalizations were 
included, which allows to apply the theory even in strongly interacting systems. 

The starting point of the derivation given here, the diffusion equation (llOli in the presence of internal 
fields, corresponds to the saddle point of the non-linear sigma model approach of 1511 . Details of the 
connection between the non-linear sigma model and the quasiclassical formalism are also given in the 
appendix. Finally we mention that in 1 53 1 a similar method was used to derive the kinetic equation in the 
presence of interaction and disorder beyond the diffusive limit. 



4 The Coulomb interaction in diffusive conductors - applications 

This section and the following are devoted to specific applications of the formalism. We will outline how 
to recover the well known Coulomb interaction correction to the linear conductivity within the present 
formalism. This will be followed by considerations on phase breaking and gauge invariance. 



4. 1 Linear conductivity 

To compute the linear conductivity we apply an external electric field E, by choosing A = — tE, <fi = 
0. Furthermore we assume that the electron distribution function has the equilibrium form, F(e, x) = 
tanh(e/2T). Equivalently, in the time domain this means F(r];x.,t) = —iT/ sinh(7rT?7). For a system 
which is homogeneous in space (after averaging over the disorder), the charge density is homogeneous, 
too, and therefore <5jBom vanishes. 

In order to calculate 5]v the screened Coulomb interaction and the diffuson are required. The dynami- 
cally screened Coulomb interaction as a function of frequency and momentum reads 

R , N 47re 2 1 -iuj + Dq 2 



q 1 u " 1 ' ^'V-iu+Dq 

The expression in the middle of this equation is valid in three dimensions. On the right hand side we 
have assumed good screening, i.e. the screening vector k, with n 2 = 87re 2 A^o, is assumed to be large. 
The perfectly screened Coulomb interaction in one and two dimensions is identical to the one in three 
dimensions, as given on the right hand side of Eq. il 141 . Next we evaluate the two diffusons D v (q, t, t') 
entering the current density. It is important to note that they appear with different time arguments ?/. In the 
second of the two diffusons in Sjy the relative time 77 is zero with the consequence that the diffuson does 
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not depend on the vector potential A, and is thus given by 

1 



-ho + Dq 2 ' 



(115) 



The convolution of the interaction with the second of the two diffusons appearing in the formula for the 
current then gives 



/ dt 2 V R (q, h - t 2 )D°{q, t 2 , t-r})= \ 1 % -t + rj), 
J 2N Dq 2 

and the expression for the current density becomes 



2e 



nT 



q 2 \smh(irTr)) 



D«{q,t- V /2,t- 377/2). 



The electric field enters via the remaining diffuson, which is given by 

D"(q,t-ry/2,t- 3^/2) = exp[-D(q- eE^r,} 

= e-^ 2 '' (1 + 2Deq ■ E77 2 + . . . ) 

After performing the momentum integration one finally arrives at 



Sjv = — 



Ae 2 D 



rid 



drj 



sinh(7rTr7) / \4irDr) 



d/2 



E. 



(116) 



(117) 



(118) 
(119) 



(120) 



in full agreement with equation ( 19 U for the correction to the conductivity. 

In the appendix it is shown how to include the electron-electron interaction in the spin triplet channel. 
The correction to the current density from the spin triplet interaction is 1 55 1 



Si 



2e 2 D 



triplet 



At 



I- A 



a\d/2 




4ttD 



drj 



■kTt) 



T Tj \smh(irTri) J J (rj - A%t) 1+d / 



dt- 



t 



E, (121) 



from which the correction to the conductivity in one, two and three dimensions are found to be 

, 2(l-Ag) 1 / 2 -2+Ag " 



<5er E EI 



-4.91 



1-3- 



^y/hD/kgT d=l 

l + 3(l + ^ln(l-4§))] ^Inih/ksTr) d = 2 , 



(122) 



1.83 



1 + 3 



2(l-Ag) 3 / 2 -2+3Ag 



d = i 



where Aq is the standard Landau parameter for interaction in the spin triplet channel. 



4.2 Charge diffusion and induced electrical potential 

One might query the fact that the correction to the conductivity due to the Coulomb interaction seems to 
be independent of the interaction strength, compare i9l\ or J 1 20b . Recalling the derivation, we note that 
this fact is related to the assumption of perfect screening. 

In this section we will present a different method |57| to calculate the screened Coulomb interaction, 
or, to be more precise, to calculate the product of the screened interaction with the diffuson. This clarifies 
further the validity of Eq. l^Kj in relation to perfect screening, and has proven to be useful also in situations 
where charging effects play a role. First notice that the diffuson describes the spreading charge cloud 
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Fig. 9 Resistivity as a function of voltage for a thin AuPd film, taken from 1581 . The resistance Ro = 11.3 MQ 
corresponds to Rq ~ 4500 fi. 



po(x, t) of a charge which is injected into the system at (xn, to). The continuity equation for this charge 
reads 



d t Po + <9x • jo(x, t) = e8(x - x )S(t - t ), 



(123) 



which is - up to the factor e - the equation for the diffuson, since j = —Dd x p. The charge density po 
induces an electric potential, and an induced charge density, which are related by the equation 



e0 ind (x,i) = y dxidtiVbarc(x, i;xi,ii)[p (xi,<i) + / o in d(xi,ii 
= J dxidtiV(x,t;xi,ti)p (xi,ti), 



(124) 
(125) 



where Vbare an d V are the bare and the screened interaction, and p; nc i is the induced charge density. The 
convolution of the screened Coulomb interaction with a diffuson, as it appears in Eqs. (Ill 31 and ( II 16> . 
is thus directly related to this induced electrical potential. In some cases the latter is very conveniently 
determined from the continuity equation for the total charge density p = po + p; n d, 



dtp + d x ■ j(x, t) = e5(x - x')5(t - t'). 
In the case of good screening this equation reduces to 

9 X • j = -aAfe d (x, t) = eS(x - x')6{t - t'), 



(126) 



(127) 



where a — 2e 2 DNo is the Drude conductivity. The induced potential as a function of frequency and 
momentum then is e</>i n d(q, ui) = 1/ (2DNoq 2 ), as found in Eq. ( II 16> following different arguments. 

4.3 Non-linear conductivity in films 



In 1979 Dolan and Osheroff [58 1 observed a logarithmic variation of the resistivity of thin metallic films 
as a function of the applied voltage; the experimental data are shown in Fig. |9] In order to explain the 
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experiment Anderson et al. | 59 1 argued that the logarithm as a function of voltage is directly related to 
the logarithm as a function of temperature (from weak localization in two dimensions) since the dissipated 
power heats the electron gas. In the case of a strong electric field, the electron temperature is of the order of 
the voltage drop on the relevant inelastic scattering length, the electron-phonon length, i.e. T ~ eEL c - p h, 
with £ c -ph = \/f Tc-ph- Assuming L -ph to be proportional to a power of the temperature, L c -ph ~ 
T~ p , electric field and temperature are related as T 1+p ~ E so that a logarithmic temperature dependence 
of the linear resistivity causes a logarithmic voltage dependence. 

Shortly after the first experiments, it was discussed whether heating is the only origin of the non-linear 
conductivity, or if an electric field - in analogy to a magnetic field - can directly destroy weak localization 
via dephasing [60 61 j. The correct answer to the second question is "no" |62 63 1, as one can easily verify 
by calculating the phase shifts of a pair of time reversed paths. In the presence of a vector potential an 
electron which propagates along a path Xi (t) or X2 (i) accumulates an additional phase 

e r v 

<Pi,2 = t / dt'Ai.a • A. (128) 

For a pair of closed time-reversed paths, X2(£) = Xi(r/ — t), and a static electric field, A(f) = — Et, the 
difference of the two phases vanishes, i.e. there is no "dephasing" from a static electric field. 

However, in a successive study of the current- voltage characteristic of gold films, Bergmann et al. | 64 1 
noted that the experimental data are not completely compatible with a pure heating model. As a possible 
explanation of the experimental findings they suggested that the Coulomb interaction correction to the 
resistivity leads to a non-ohmic contribution, namely 

e 2 

6R EE1 /R 2 = -^—\n[(k B T) 2 + aHD(eE) 2 / k B T], (129) 

where a is a factor of the order one. 

Indeed in formulating the phase shift argument for the interaction contribution, a sensitivity to a static 
electric field cannot be excluded: The interaction correction to the conductivity is related to the propagation 
of a particle and a hole along closed paths. In the absence of a vector potential the phases of particle and 
hole cancel, whereas in the presence of a vector potential the accumulated phase difference is 

Sp=j di'xi A - - dt'x 2 ■ A. (130) 



KJti-r, h ■ . 

As is shown in Fig.[H]the relevant paths obey the relations Xi(t) = x%(t) for < t < t\ and X\(t — rj) = 
X2(t) for t\ < i < rj, which allows the phase difference to be written as 

e r° 

5<p = - ttt'ii ■ [A(i') - A(t' + r])]. (131) 

" Jt-l-T) 

For the particular case of a static electric field described by A = — Ei, this phase shift becomes dip = 
|?y(x2 — xi) ■ E. This suggests that the interaction correction should be sensitive to a static electric field, 
leading to a non-linear conductivity. The quantitative calculation in fact results in 1271 

e 2 / hD(eE) 2 \ 
6R EEI /R 2 a --— lnk B T + 1.62 V ' , (132) 
2-K z n \ (TTKB 1 ) J 

verifying the non-ohmic behavior of the resistivity with the characteristic electric field scale as suggested 
in |64|. In order to derive Eq. Jl 32b . a thermal distribution function with electron temperature T has been 
assumed. Following the steps in our calculation of the linear conductivity in section |4~TI Eq. (II 32i is ob- 
tained from il 17> in two dimensions under the condition hD(eE) 2 <C (nksT) 3 . Note that the functional 
form ( f!29i of the non-ohmic resistivity for strong electric fields cannot be confirmed theoretically. We 
remark, finally, that further theoretical results have been given in 1271 1551 . related to strong electric fields, 
time dependent fields, one to three dimensions, and magnetic field effects. 
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4.4 Gauge invariance 

In the previous section we demonstrated that the non-linear electric field effect in the resistivity can be 
understood in terms of dephasing by calculating the phase shifts of the relevant classical paths in the 
presence of a time dependent vector potential. One may argue against this interpretation by noting that 
a vector potential can be gauged away in such a way that the static electric field is described by a static 
scalar potential, E(x) = — <9 x 0(x). A static scalar potential does no longer affect the diffuson, so that the 
argument of the phase difference along two classical paths can no longer be used. In order to clarify the 
situation, we now verify explicitly the gauge invariance of the expression for the current density. First one 
notices that <5jBom = —Dd x 5p is obviously gauge invariant. For 5jv, on the other hand, an explicit check 
is necessary. Given the gauge transformation 

A - A + d xX (133) 
4> -> <t>-dtx (134) 
the distribution function and the diffuson transform according to 

Ffa;x,f) -» F(t/;x,t)exp{-ie[x(x,t+^)-x(x,t-|)]} (135) 
L> r '(x,i;x',t') -> D*'(x,t;x',t / ) 

x exp{-ie[ X (x, t + |) - *(x, t - |)]} 

x exp{ie[ X (x', t' + |) - X (x', t' - |)]}. (136) 

By applying the above transformation to 8]v as given in Eq. J 1 1 3I >. one can easily verify that the function 
x(x, t) drops out, so that the expression is manifestly gauge invariant. 

For the specific example of a static electric field with A = — Ei we choose x( x 7^) = tE • x. After 
the gauge transformation the electric field appears in the distribution function, tanh(e/2T) — * tanh[(e — 
/z x )/2T], /i x = eE • x, but not in the diffuson. We conclude that although the contribution of the Coulomb 
interaction to the current density is gauge invariant, the interpretation of the non-linear electric field effects 
depends on the actual choice of the potentials A and <fr. With E = — d t A, cf> — we would interpret the 
non-linear conductivity as due to a phase effect. With E = —d x (j>, A = the origin of the non-linear 
conductivity can be attributed to the varying local chemical potential which is felt by the traveling particle, 
and the virtual hole created and absorbed in the intermediate states. 



4.5 Electron dephasing 

The dephasing time sets the scale over which an electron propagates without loosing phase coherence, 
and determines the amplitude of the quantum interference corrections to the conductivity. Since elastic 
scattering does not destroy the electronic phase coherence, inelastic scattering must be hold responsible. 
Dephasing has been studied in much detail for weak localization 1 65 66 67 68 1 (for a recent experimental 
review see [69 1), and also in the context of universal conductance fluctuations l70ll7Tl . The amplitude of 
the interaction correction to the conductivity, on the other hand, is set by the thermal time tt = h/kBT. 
In most cases the dephasing time is much longer than the thermal time. In some experiments, however, an 
extraordinarily strong phase breaking has been reported: Recently a low temperature saturation of t v in 
gold wires 1 72 73 1 has attracted much attention 1 23 , 24 , 25 28 74 75 76 29 1. Furthermore, in a number of 
cuprates, the dephasing rate decreases only slowly with decreasing temperature 1771 1781 l79l . For example 
in Bi2Sr2CuC>6, 1771 . the dephasing rate varies as 1/t v ~ T 1 / 3 , with t v much shorter than tt- In such a 
case phase breaking may become relevant also in the interaction contribution to the conductivity. 

Unfortunately, in the experiments cited above, it is not clear which microscopic mechanism is responsi- 
ble for the strong phase breaking. Nevertheless we believe it is important to answer the following questions: 
(i) Is phase breaking relevant in the interaction contribution to the conductivity? (ii) If yes, is the phase 
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breaking rate which is relevant in the interaction contribution to the conductivity the same as for weak 
localization? While Castellani et al. | 80 1 came to the conclusion that the answer is "yes" to both questions, 
Raimondi et al. 1271 who reexamined the problem, confirmed the suggestion that dephasing is relevant 
also in the interaction contribution to the conductivity, but found a different dephasing time for the latter 
compared to weak localization. After some historical remarks we will give a semi-quantitative summary 
of the analysis carried out in |27|. 

While it is clear that inelastic scattering contributes to dephasing, the precise way this happens is less 
obvious. It was first noted by Schmid 1 8 1 1 that the inelastic quasi-particle scattering time is enhanced in the 
presence of disorder. The dephasing time was initially assumed to be identical to the inelastic quasi-particle 
scattering time |82|, thus the inverse dephasing time was assumed to be oc TlnT in two dimensions. 
Some time later, Fukuyama and Abrahams 1 66 1 reexamined the problem in terms of standard diagrams and 
calculated the "mass" that develops in the particle-particle propagator (the cooperon) in the presence of 
Coulomb interaction. They also found an inverse time proportional to T In T. Furthermore a mass term also 
appears in the particle-hole propagator (the diffuson) 1 80 1 . The mass in the particle-hole propagator turned 
out to coincide with that in the particle-particle propagator. However, calculating a mass is not sufficient 
in order to draw final conclusions; Altshuler, Aronov, and Khmelnitskii |65|, for example, calculated the 
dephasing time via a path-integral approach and predicted an inverse dephasing time proportional to T, in 
contrast to |82|. This result - see also |67| - has been confirmed in many experiments. The approach of 
Altshuler, Aronov, and Khmelnitskii has been extended later in order to determine the relevant dephasings 
times for the universal conductance fluctuations 1 7 1 1 , and in the particle-hole channel 1 27 1 . 

The general idea is that the phase of an electron which is propagating along a path xi^t) is shifted by 
yi,2 = e J dt xi 2(t') • A(t'). In the presence of fluctuating internal fields this results in phase fluctuations 
which finally destroys phase coherence. In thermal equilibrium the internal electric field fluctuations are 
given by 



where V (q, ui), V R (q, oj) are the Keldysh and the retarded Coulomb interaction, and low frequencies 
(w <C T) have been assumed. The assumption of low frequencies is essential here since only in this 
frequency regime the classical field fluctuations, i.e. the Keldysh component of the interaction, dominates 
whereas at higher frequencies the nontrivial structure of the interaction in Keldysh space has to be taken 
into account. It has been shown by Fukuyama and Abrahams |66| and Castellani et al. |47| that the high 
frequency contributions of the Coulomb interaction renormalize parameters like the diffusion constant but 
do not contribute to dephasing. It is therefore reasonable to neglect these high frequency contributions in 
the following. The particle-particle propagator (cooperon) and the particle-hole propagator (diffuson) in 
the presence of the fluctuating field can then be written as a path-integral: 



e 2 {E l E j ) 



(137) 



-q i q j — ImV fl (q,u;), 



2T 



(138) 




(139) 




(140) 



with 




(141) 



(142) 
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In both cases the Si term describes the additional phase due to the vector potential A. Since we have 
to consider the phase difference of pairs of classical paths, the vector potential appears twice in Si. In 
the case of the cooperon the two paths are the time reversed of the other, i.e., to a velocity x on path 
"1" corresponds a velocity — x on path "2" and in the expression for the phase difference the sum of two 
vector potentials appears. In the case of the diffuson the two paths are traversed in the same direction, but 
at different times. Therefore the minus sign between the two vector potentials remains, and in particular 
when the two paths are traversed at the same time (77 = 0) the phase difference vanishes. By averaging the 
phase factor exp(— Si) over the electric field fluctuations, using the relation 



(<<xi)(i-V)} = <'M> I -\{{5<p?) 



exp(-S), 



(143) 



one finds 



i/>/>/;j??[-^(-) 



exp{iq • [xi(ii) - xi(t 2 )]} 



cos uo 



h-t 2 



ti+t 2 



(144) 



in the case of the cooperon, and 



S?- h - = fdtjdtj ^.^—[-hnV^)] 

JO JO J-T Z7r u 

x exp{iq- [xi(t 1 )-x 1 (t 2 )]}e- ia 'C* 1 -* 2 ) [1 - cos(u;t?)] 



(145) 



in the case of the diffuson; in both S pp ' and S p the velocities x(ii, 2 ) were removed by a partial inte- 
gration, and the relation between the two relevant paths xi (t) and x 2 (t) was exploited. The expressions 
for S p p ' and S p appear to be similar, except for the time dependent factors in the second line; this 
leads to different dephasing times in the particle-particle and particle-hole channels. In two dimensions, 
for example, the final results are 



S p. P . 



h/e 2 



h/e 2 



Rr 



(146) 



and 



s ph - 



(T^-^HDk 2 ^ 2 ); 



Dk 2 (h/e 



T 



Ru 



(147) 



where k is the inverse screening length, and Rq is the sheet resistance. The dephasing time is here deter- 
mined from the condition S(i~i = t v ) = 1. In the particle-particle channel the standard result of Altshuler, 
Aronov and Khmelnitskii 1651 1671 168 1 is confirmed. The dephasing time in the particle-hole channel is 
finite, different from the one in the particle-particle channel, and is identical to the inelastic scattering 
rate in the two-particle propagators (diffuson or cooperon), as first calculated by Fukuyama and Abrahams 
|66|IH4|. 

The low frequency electric field fluctuations considered here, cannot explain the strong phase breaking 
observed in 1721 1771 . However we expect similar results for other phase breaking mechanisms: We sug- 
gest that a mechanism leading to strong dephasing in the particle-particle channel will also cause strong 
dephasing in the particle-hole channel. In particular when t v becomes comparable to or shorter than the 
thermal time tt> the relevant time scale which sets the amplitude of the Coulomb interaction contribution 
to the conductivity will be t v instead of tt- This is consistent with the experiments: In the gold wires of 
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1 72 1 the dephasing rate saturated below T ~ 1 K to values of the order of H/t^ ~ 1-10 mK. In the samples 
with the strongest phase breaking a saturation of the interaction correction to the conductivity has been 
observed below T ~ 100 mK (73). In B^S^CuOe, a compound with a single CuC>2 plane per unit cell, 
the in-plane zero field resistivity increases as In T below ~18 K, consistent with quantum interference ef- 
fects in two dimensions. The shape of the orbital magneto-resistance is well fitted by the weak localization 
expression [77 1, but with an unexpected large dephasing rate which varies as T 1 / 3 . The spin component of 
the magneto-resistance varies at low fields as [R(B) — R(0)]/R(Q) ~ (B/B s ) 2 . From the standard theory 
for the Coulomb interaction in the spin triplet channel 183111 51 one would expect such a quadratic magnetic 
field dependence with a magnetic field scale which is linear in the temperature, B s ~ 1/tt. Experimen- 
tally the magnetic field scale varies as T° 4 , which is close to the temperature variation of the dephasing 
rate. This suggests that the origin of the spin component of the magneto-resistance might indeed be the 
Coulomb interaction contribution to the conductance, but in the presence of an up to now unidentified 
phase breaking mechanism. 



5 The Coulomb interaction in nanostructures 

Quasiclassical methods - as pioneered by Nazarov - allow description of interaction effects in a unified 
way. For instance, in the following subsection we will show that Coulomb blockade physics in small 
tunnel junctions and disorder-induced zero-bias anomalies in diffusive conductors are closely related. More 
generally the question arises how the Coulomb interaction affects the transport in small systems with 
diffusive electron motion, where interfaces provide an additional source of scattering. Typical situations 
will be examined in detail in the subsequent subsections. 



5.1 Coulomb interaction and interface conductance 



The Coulomb interaction plays an important role not only in diffusive transport, but also for transport 
across interfaces. For example, the Coulomb interaction reduces the density of states near the Fermi energy, 
1201 12 1 1 1 1031 1461 1471 1851 1861 15 1 1 which leads to a suppression of the transmission of tunnel junctions, and to 
a non-linear current-voltage characteristic at low bias. On the other hand Coulomb blocking of tunneling 
may also be due to interface charging effects |87 88 89 90 91 1. These two classes of phenomena were 
described in a unified way for the first time by Nazarov |92 93 1, and later by various authors |85 95 96 
l97ll98ll99lfT00l. 

In this subsection we will show how to construct the theory within the quasiclassical scheme. Whereas 
in the above cited references it is assumed that the electronic system on both sides of the interface is in 
thermal equilibrium, there is no such assumption in the approach we present here. The idea is to start 
with the boundary conditions for the quasiclassical Green's function, derive an expression for the current 
through the interface, and at the end average over the fluctuations of the internal electric field. 

Zaitsev's boundary condition 1 10l1 11021 . which generalizes the matching condition (1571 to situations 
where g R ' A ^ ±1, reads 



T 

R - Raa + — (s(x_) - s(x + )) 2 



T 

= j[s( x +M( x - 



(148) 



with 



<Kp) + <?(p) 

<?(p) - ff(p) 



(149) 
(150) 



The antisymmetric combination of Green's functions a is continuous across the interface. In the absence 
of electronic interactions (g R ' A = ±1), we recover Eq. (I57> . In the interacting case, where g R,A ^ ±1, 
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we restrict our considerations to the limit of low interface transparency, T/R -C 1 . Then the structure of 
the boundary conditions becomes simple, 

a K = ^[s(x+),s(x_)] K , (151) 

and the current density across an interface which connects two dirty pieces of metal is given by 

1 f d 

n-j = enN T- J ^([5 s (x + ),.9 s (x_)] K ) 0i „t (152) 



1 



T = v F - / d(cos(9)cos6>T(cos6>). (153) 
2 Jo 

In the absence of interactions the tunneling conductivity of the interface is given by ot = 2e 2 NoT. The 
remaining technical difficulty in order to evaluate the current through an interface in the presence of in- 
teractions is the average over the internal field fluctuations. Solutions are known only in some special 
situations; one of these is when the two sides of the interface are uncorrelated such that the average over 
the internal field fluctuations factorizes. The current can then be written in a familiar form in terms of the 
local density of states N(e, x) as 

n-j = /de[F(e,x_)-F(e,x+)]iV(e,x + )iV(e,x_) (154) 

iV(e,x) = ±N Re(g R ' A (e,x))^. (155) 

In general, however, such correlations do exist due to the interaction between the electrons on opposite 
sides of the interface. 

Another situation where Eq. dl52> can be evaluated explicitly is when the gradients in the kinetic equa- 
tion on both sides of the interface are negligible, i.e. when the diffusive motion of the charge carriers is 
negligible. In this limit standard Coulomb blockade theory (P(E) theory, 1571 1 is recovered as we demon- 
strate here. The question under which circumstances the gradient terms in the kinetic equation can be 
neglected will not be addressed here: For a discussion of this issue see the appendix of |57|, and also sub- 
section Neglecting the gradient term, the kinetic equation in the presence of the Coulomb interaction 
can be solved explicitly, and the Green's function reads 

3 tlt2 (x) = e^^so.WsMe-^'*), (156) 

with dt<f> = e^int, and go is the Green's function in the absence of the interaction. Since we assume 
Gaussian fluctuations for the internal field 0i nt , it is clear that the average over </>; nt can be carried out for 
both the Green's functions and the current density, Eq. dl52t . The final result for the current density can 
be written in the form 

n-j = e2IW J de J dwP(w) 

x {/(e, x+)[l - f(e - u, x_)] - /(e, x_)[l - /(e - w, x+)]} (157) 

withF(e,x) = 1— 2/(e,x). The factors /(e,x+)[l — /(e— w,x_)] describe the tunneling of acharge from 
an occupied state on the "+"-side of the junction to an unoccupied state on the "—"-side. In the absence of 
interactions the function P(lo) is a delta function centered at zero, and the usual expression for the tunnel 
cunent is obtained. In the presence of interactions a tunneling electron may exchange energy with the 
environment, and P(u>) may be interpreted as the probability to transfer the energy lu to the environment 
1571 . The function P(lo) is given by 

P(to) = -!- / dte iuJt e J{t) (158) 

J(t) = ((v'i(«)-V'i(0)M(0))^ rt +<V'i(«)^(0)>^ < ,-(^(t)v'i(0)>^ > (159) 
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where tp 12 is the phase difference between the left and right leads, <p\£ = <pi,2 (x_) — ^1,2 (x+); the index 
1, 2 refers to the Keldysh space. In thermal equilibrium and expressing the phase/voltage fluctuations over 
the junction in terms of an impedance Z one has 

dte^{ Vl (t)M0)h iat = \—\ (160) 
the function J(t) is determined as 

J{t) = J^^p{ [cosiujt) ' 1]coth {^)' is[n{ujt) }- (161) 

As an example, to mimic the electrodynamic environment of a tunnel junction, assume the impedance Z to 
be Z = l/(— iuiC + R^ 1 ), i.e. a capacitance in parallel with a resistor. For a high impedance, R^S> h/e 2 , 
when the charge relaxation process becomes slow, one obtains 

J(t) = —iE c t — E c Tt 2 , (162) 

where E c = e 2 /2C is the charging energy. In the zero temperature limit P(u>) reduces to a delta function, 

P(u) =S(uj-E c ), (163) 

such that a minimum voltage difference E c /e is needed for an electron to tunnel through the junction. The 
current- voltage characteristic is given by 

I(V) = G T {V - E c /e)G(eV - E c ), G T = 2e 2 N TA. (164) 

For more applications of the standard Coulomb blockade theory we refer to the literature, e.g. II91I . 

We close this subsection by highlighting how the zero-bias Altshuler-Aronov anomalies arise. These 
are present in situations where the gradient in the kinetic equation cannot be neglected. Since under this 
condition an exact expression for the Green's function near the junction is not known, we restrict ourselves 
to a pertubative expansion in the internal field 4>- m t, with the result 



n-dj = - e ^V ri||i[F(e,x_)-F(6,x + )] 

x Re( ( 5. g K (x + )[<5/(x + ) - ^(l)])^ 
+Re(6g K (x-)[Sg z (x-) - <5 9 z (x + )]) ,„ t 



(165) 



where 5g K and Sg z are the correction to the Green's functions due to the fluctuating field. The relevant 
correlation function has been given in Eq. ITlOi . and we find 

nSj = J de[F(e,x_)-F(e,x + )]Im J^J dx x ^(e - w,xi) 

x {p (x+,w;xi) [0 ind (xi,w;x + ) - ind (xi, w; x_)] 

+ j o (x-,cj;x :L ) [</> ind (xi, w;x_) - in d(xi, u; x+)]} , (166) 

where ctt is the tunneling conductivity in the absence of interactions, po(x, x') and <^i n d(x, x') are the 
charge density and the induced electrical potential at x due to a charge which is generated at x', compare 
subsection l4.2l When interactions are taken into account only on one side of the interface, the correction 
to the current is controlled by the density of states. For example in three dimensions the correction to the 
density of states is proportional to yfe, giving rise to a \/V of the tunneling conductivity as a function of 
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Fig. 10 The geometries considered in subsections 15. 21 — [5.41 Section l5~2l is devoted to wires which are much longer 
than the thermal diffusion length Lt- In Sec. l5.3l i/T is of the order of the system size; both in Sec. l5.2l and Sec, l5.3l we 
assume ideal leads. In section l5"4l a diffusive point contact is considered, where the system size is much smaller than 
Lt such that the phase coherent volume extends into the leads, which then have to be taken into account explicitly. 

voltage 1131 . However, in the more general situation, as pointed out by Nazarov 1921 1931 . correlations in 
both leads and between the leads have to be taken into account. 

For further experiments where the perturbative zero bias anomalies in the tunnel conductivity have been 
observed, see for example 1 96 97 98 1. We will work out the Coulomb correction ( I166> to the interface 
current for a specific geometry in section 1531 



In this and the following two subsections we will discuss the contribution of the Coulomb interaction to 
the conductivity in small systems. We would like to mention that size effects in the quantum corrections to 
the linear conductivity were studied in detail in the mid-eighties 1 1031 1 1041 1 1051 . In contrast, we will focus 
here on the case of a large applied voltage. We will successively decrease the size of the systems under 
consideration. This subsection is devoted to wires which may be shorter than the inelastic scattering length, 
but which are still longer than the thermal diffusion length Lt — \J hD j ksT . This is then followed by 
subsection !5. 31 where we consider system sizes of the order of Lt- There it will be essential to consider 
the discreteness of the diffusive modes, and we will also allow for resistive interfaces between the wire and 
the leads. Finally, in subsection !5.4l we will consider the situation of a small, fully coherent piece of metal 
beeing connected to two large reservoirs. 

We start with a long, thin wire, for which the Coulomb correction to the resistivity, including the leading 
non-linear voltage dependence, in analogy to Eq. dl32t for a thin film, is given by 1 55 1 



where V = EL is the applied voltage. Again, this result has been obtained under the assumption of a 
thermal distribution function with a constant temperature T. Whereas this is reasonable for macroscopic 
samples, it fails, however, in samples which are shorter than the electron-phonon scattering length, in 
which case the Coulomb interaction correction to the conductivity far from equilibrium |56 106 , 54 1 has 
to be considered. 

For an evaluation of the current-voltage characteristic, the diffuson and the distribution function in the 
wire are required. The diffuson is found by solving the differential equation dl08t with the condition that 
the derivative normal to an insulating boundary (i.b.) vanishes, i.e. 



5.2 Non-linear conductivity in wires 




(167) 



(n- V x )D(x,x') 



x G i.b. 
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while the diffuson itself vanishes at a metallic boundary (m.b.), 

^ x ' x ')| xem . b .=0. (169) 

The first condition means that no current flows through the interface, the second condition corresponds to 
the assumption that an electron arriving at the metallic boundary escapes into the leads with zero probability 
to come back into the wire. Furthermore it is assumed that the left and right leads of the wire are in thermal 
equilibrium, 

^M)L su .,,=-(«). ("0, 

By solving the kinetic equation it is found that the distribution function depends on the various relaxation 
mechanisms governing the collision integral 1 107 108l ll09lll 101 1241 . and we distinguish three regimes: 

a) When the system is much longer than the electron-phonon scattering length, L ^> L _ p h, the distri- 
bution function acquires the equilibrium form with a local chemical potential and temperature, 

where x = ... Lis the distance from the left lead. The local electron temperature T e (x) may be de- 
termined from an energy balance argument, assuming that the dissipated power equals the gradient of 
the heat flow. In the limit considered here the heat flow is dominated by the phonons. For a stationary 
temperature T e (V) the Joule heating power P ln = erE 2 equals the power which is transferred into 
the phonon system, P ou t- For weak heating one has P out = cyAT/r c _ph, where cy is the electron 
specific heat, AT is the difference between the electron and phonon temperatures, and T c _ p h is the 
relevant energy relaxation rate. For strong heating, on the other hand, the effective electron tempera- 
ture is of the order of the voltage drop over a phonon length, T e ~ |e| VLc^p^/ L. If the "hot" phonons 
escape ballistically into the substrate, the temperature in the bulk of the wire does not depend on the 
position x but is voltage dependent, T e — T e (V). By neglecting the region near the leads, where the 
temperature rises from Xi oa( j to T e (V), Eq. dl67> is recovered for the voltage dependent resistivity. 

b) When £ c -ph 3> L ^> L m one still expects a distribution function near local equilibrium, due to 
electron-electron scattering. The local temperature is determined from the relation aE 2 =— V[/cVT e (x)], 
where k is the thermal conductivity. Using the Wiedemann-Franz law, k — (7r 2 /3) T e cr(fcs/e) 2 , the 
temperature profile in the wire is determined as 

Ti(x)=T^+^(^pj x(L-x). (172) 

c) In the absence of inelastic scattering, realized when L ln 3> L, the distribution function is a linear 
superposition of the distribution function of the leads, 

F(e, x) = [{L - x)F{e, 0) + xF(e, L)] /L. (173) 

The resistance as a function of voltage in regime a) is given in Eq. J167I . This equation applies when the 
voltage drop over a thermal diffusion length is smaller than the temperature, eVLx / L < T. Since we 
assume that £ c -ph > L the electron temperature as a function of voltage rises so fast that this condition 
always holds. We compare now the relative importance of both heating and non-heating effects for the non- 
linear resistivity. From the linear conductivity and the increase of temperature due to the applied voltage 
as discussed above, we find at low voltage 

SR^JR ~-4.91^ — T2(1/Tc _ phr ( 174 ) 
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Fig. 11 Interaction correction to the conductance SIeei/V for a mesoscopic wire as a function of voltage, taken 
from |54|. SIeei/V is plotted in units of (e 2 /K)Lt / L. The full line corresponds to the non-equilibrium distribution 
function 11731 . The line with long dashes corresponds to the local equilibrium distribution function 11711 with x- 
dependent temperature. The line with short dashes (L/Lt = 5) is the non-linear conductivity due to the heating 
contribution only, Eq. J 1781 . 



which has to be compared with the corresponding change of the resistance due to non-heating effects: 

2 e 2 L T D(eV/L) 2 

dit n on-heating/-R ~ — 0.21— r— — r . (175) 

TT L 1 A 

One observes that the heating contribution to the non-linear resistivity is larger by a factor Tr e _ph. Thus 
the non-heating effects are difficult to observe experimentally, since usually l/r c _ p h -C T. 

We now consider the other two limits, b) and c). We found that in both limits the current can be written 
as 1 54 1 

M"eei(V,T) = ^~S {eV/k B T) V, (176) 
n L 

where the function / depends on the distribution function and on the length of the wire; T is the temperature 
in the leads. Numerical results are shown in Fig. ^3 Notice that f (eV / ksT) is proportional to SIeei/V, 
so that /(eV/fcsT) also represents the voltage dependent conductance in units of {e 2 /K)(Lt / L). For low 
voltage and large system size L ^S> Lt, the result 1 103 1 



SIeei/V ^ -^-^f 
n L 



3C(3/2) 5 L T 5 % /2C(5/2) L\ 
(2tt) 3 / 2 ~ 6 T + ^ T 2 



3C(3/2) 

- 0-498 (177) 



is approached. With decreasing length the correction to the current decreases, since electrons escape more 
quickly from the wire into the leads. The full lines show the voltage dependent conductance for case c); 
the long-dashed line corresponds to case b). The short-dashed line (L /Lt — 5) is obtained within a simple 
approximation: Instead of evaluating the full expression for (5]eei we take the linear conductivity as a 
function of temperature, and average over the temperature profile, 

^heating = -j- / dxSa(T(x)). (178) 
L Jo 

Important results are the following: in both cases b) and c) the conductance scales with voltage over 
temperature. In case b) (hot electrons) the main effect is simple heating, i.e. the non-ohmic effects are small. 
In case c) (far from equilibrium) the current-voltage characteristic is quantitatively different from the hot 
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Fig. 12 Temperature dependence of the Coulomb interaction correction to the conductance of a mesoscopic wire in 
arbitrary units. The full line corresponds to the double-step distribution function, the dashed line corresponds to hot 
electrons with temperature profile, as explained in the text; T is the electron temperature in the reservoirs. 

electron regime. The temperature dependence of the Coulomb interaction contribution to the conductance 
is shown in Fig.[21in a double logarithmic plot. The curves shown are obtained from the same numerical 
data for /(eV^/fcsT) as shown in Fig.[S\(L / Lt = 200). At high temperature the Coulomb correction to 
the conductance follows 1 / \/T, which is seen as a linear behavior in the double logarithmic plot. When 
the temperature in the leads becomes lower than the voltage, the conductance saturates. In the absence of 
inelastic scattering, c), this low temperature saturation appears at a higher temperature than in b). 

Finite size effects in the linear conductivity of thin wires of length L ~ have been studied e.g. 
by Masden and Giordano |104| who found a qualitative (although not quantitative) agreement with the 
theoretical predictions. The distribution function in mesoscopic wires out of thermal equilibrium was 
measured by Pothier et al. II 1 II . For short wires and low temperature (L « 1.5 /im, T w 25 mK), the 
double step like distribution function was observed. Unfortunately, as far as we know, there is no detailed 
investigation of the temperature and voltage dependence of the conductivity in this experiment. 

5.3 Short wire with interfaces 

In subsection l2.4l we demonstrated how to obtain the classical resistance of a system which is composed 
of diffusive pieces and resistive interfaces in the framework of the quasiclassical description. In this sub- 
section, following closely 111 121 . we include the Coulomb interaction, and in particular we will concentrate 
on structures of a size comparable to the thermal diffusion length Lt- In subsection 15. II we described a 
general formalism to obtain the quantum correction to the current through an interface; here we apply the 
formalism to a specific geometry where a diffusive piece of metal is connected via resistive interfaces to a 
left (x = 0) and right (x = L) reservoir. 

We divide the quantum correction to the current into two contributions: The first has the meaning of a 
correction to the conductance, 51^ = V5G. The second can be interpreted as the redistribution of the 
voltage across the interfaces and the wire, 51^ = G5V. More precisely we denote by 51^ those terms 
which arise from the correction to the distribution function like, for example, 



we will give the explicit expression for SI^J below. Both terms are necessary in order to ensure current 
conservation. In principle the variation of the distribution function, 5F(e; x), can be determined by solving 
the kinetic equation in the presence of interactions. However by exploiting current conservation it appears 




(179) 
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that this is not necessary. For the time independent problem which we consider here, current conservation 
means that the current does not depend on the position, i.e. Ii = I W irc = Ir, where the subscripts I and r 
indicate the left and right interface. By fixing the voltage drop over the whole system to V, such that the 
sum 5Vi + 5V r + SVwiTc is zero, it is possible to eliminate 51^ from the equations to obtain the correction 
to the current as 

dlvm = p p 5 • (180) 

Hi + H r + il w i re 

To proceed further we need the explicit form of 8IiJ- Ie . For an interface which is attached to an ideal 
lead on one side, the quantum correction to the current is controlled by the correction to the density of 
states on the other side, 

51^ = -eATi J deSN(e,0)[F(e;0) -F l (e)} (181) 

SlP = -eAT r J de5N(e,L)[F r (e)-F(e;L)]. (182) 
Quite generally we obtain the density of states correction as 

/J 
■^S(x,x) (183) 

S(x,x) = — Im J dxi 

xF(e - uj; Xx)p Q (x, uj; x 1 )^ ind (a; 1 , u>\ x), (184) 

where po(x,u>; x%) describes the spreading of a charge injected into the system at xi, and 0j n d is the 
induced electrical potential as we discussed in subsection l4.2l By using the same notation for the correction 
to the current in a diffusive wire we find 

= eDN AjdeJ^-d x [F(e;x)S(x,x)} (185) 

~2eDN A J de J ^F{e; x)d Xl S(x, Xl )\ Xl=x . 

In general 5I^ Ie (x) depends on the position x. Equation (1 1 801 . however, is constructed in such a way that 
the spatial average SI^ = L^ 1 J dxdl^ (x) has to be inserted. 

We have already discussed the distribution function F including the relevant boundary conditions in 
subsection l2.4l Inside a diffusive wire the charge density po satisfies the diffusion equation, Eq. dl23> . At 
the boundaries to the left and right reservoir, one may derive the matching conditions 

Dd x p (uj,x;x')\ x =o = Tip o (uj,0;x') 

Dd x p (u>,x;x')\ x = L = -r r p (u,L;x'), (186) 



to be compared with the boundary conditions for the distribution function, Eqs. (1621 and i63\ . A careful 
analysis is also required for the induced field 4>- m d(ui, x; x 1 ). We assume good metallic screening (compare 
subsection 14. 21 so that an injected charge is almost instantly screened, and therefore the wire will be 
electrically neutral with the possible exception of a thin surface layer. In this case inside the wire, one has 

— ad x ^i a d((^, x; x') — eS(x — x'), (187) 

where a — 2e 2 DNoA is the conductivity. 

From now on we concentrate on a short resistive wire with interfaces assuming temperatures of the order 
of and lower than the Thouless energy, HD/L 2 . We start by expanding the charge density po in diffusive 
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Fig. 13 Contribution of the electron-electron interaction to the temperature dependence of the conductance of a wire 
with length L ~ Lt for different values of the bias voltage, from 1 1 12 1. The results are obtained under the assumption 
that the charge diffusion can be described with just a single diffusive mode. dl/dV is in units of e 2 /h and must be 
multiplied with the non-universal number A defined in the Eq. 11901 - 70 is the energy of the lowest diffusive mode; in 
the absence of interface barriers 70 = iv 2 hD/L 2 , for strongly resistive interfaces 70 — > 0. 



modes, 



po(u),x;x') = 



fn{x)fn(x') 



(188) 



-lui + 7„ 

where the (normalized) functions /„ (x) are obtained from the eigenvalue equation 

-Ddlf n {x)= ln f n {x). (189) 

In the zero-dimensional limit we approximate the sum in Eq. i ll 881 by retaining only the eigenmode with 
the lowest energy, i.e. p u (x,x') — ► efo(x)fo(x')/(—iui + 70). This approximation is justified when 
the energy scales related to the temperature and to the voltage remain below the energy of the second 
lowest diffusive mode. As a consequence of Eq. dl 871 the field (f>i n d(uj, x; x') is frequency independent 
and one observes that the frequency dependent factors in all the contributions to the current are identical, 
6I(T, V) ~ Fl_ u Fl/(-\w + 70). The explicit result reads 



SL 



EEI 



-A— / drye' 70 '' 
2Wo 



sinh(7rT7j) 



sin(eyr7), 



(190) 



where only the dimensionless number A and the quantity 70 depend on the details of the system under 
consideration. A short time cut-off has to be introduced in the integral in order to avoid a logarithmic 
divergence. An explicit form can be obtained by refining the approximations made. However we do not 
specify the cut-off here since it only leads to a temperature independent correction to the linear conduc- 
tance. 

Let us first discuss the temperature and voltage dependence of SIeei, then we determe A and 70 explic- 
itly in the two limits of perfectly transparent interfaces and for interfaces with low transparency. FigurefOl 
shows 81/ dV I (Ae 2 /2tt) as a function of temperature, a constant being subtracted. At high temperature a 
logarithmic behavior is observed, 



^ = G(T) = G(T ) 



Aln(T/T ), 

2irn 



(191) 



which saturates below T sat ~ max(7o, eV). Figure^Jshows the voltage dependence of the conductance; 
here the linear conductance has been subtracted. For 70 <C T the conductance scales with voltage over 
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Fig. 14 Voltage dependence of the conductance of a wire of length L ~ Lt- In all curves the linear conductance has 
been subtracted. The figure is taken from 11121 . 



temperature, while when 70 is large the relevant scale for conductance variations is 70. Notice that the 
logT behavior and the current-voltage characteristic are universal for the zero-dimensional systems: The 
results are not sensitive of whether the system is dominated by interfaces or diffusion, nor to the shape of 
the sample, nor to the precise value of the conductance. The microscopic details are only relevant via the 
number A and the escape rate 70. 

How large are the amplitude A and the energy of the lowest diffusion mode 70? In the case of two well 
transmitting interfaces, G W i ro <C G; , G r , the eigenfunction of Eq. d!89i with the lowest eigenvalue is 

fo(x) = yjsin(7rz/L),7o = n 2 D/L 2 . (192) 

The distribution function and the potential ^; n( j are determined as 

F(x) = [{L^x)F l +xF r ]/L (193) 

± 1 i\ e f (L~x')x/L 2 x<x' , 1n .. 

= G— { \l-x)x',L 2 x> X > (l94) 

and the amplitude of the correction to the current is found to be A — 64/7T 4 — 4/ir 2 w 0.25. In the 
opposite limit, G w i re 3> G; , G r the eigenvalue equation dl 89t may be solved perturbatively in the barrier 
transparency and one obtains 



fo(x) = l/Vi,7o = (rj+r r )/L (195) 

F(x) = (r i F i + r rJ p r )/(r i + r r ) (196) 

<f> ind (u,x;x') = e/id+Gr), (197) 

which leads to A = 2r J r r /(r J + T r ) 2 . 

It is useful at this point to briefly discuss how the charging effects modify the above results. We assume 
that the accumulated charge is proportional to the field as p(x) = po + Pmd = C4>i n a(x), where G is a 
capacitance per unit length. From the continuity equation for the total charge density p(x), one obtains a 
diffusion equation for the field 

— kjG0md(w, x; x) — adl<j)i n d(ui, x; x') — eS(x ~ x'). (198) 

For highly transmitting interfaces the solution is 

^ d ( W ') = 2eE Sm( 7;/fl S ; n( ^ /L) . (199) 

i! — ' —1U){CL) + (7rn)^G w ire 
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and the correction to the current is modified according to 



£>/eei = -— / drje 

*K Jo 



ttT 2 



svn(eVri) 



sinh(7rT?7) 

x ^2 A ni l ~ exp[-r ? (7m) 2 /i?(CL)]}. (200) 

n 

The numbers A n depend on the wave-function of the diffusive mode n, and R = G~l TC . Notice that 
although the free charge density po(x) is zero-dimensional many diffusive modes have to be taken into 
account in the field <^>i n< j (x) . For a small capacitance C the charging simply cures the short time divergence 
in the integral in Eq. ( I190i . For a larger capacitance the /-^-characteristic is no longer a universal function. 

For the system with poorly transmitting interfaces, the field in Eq. dl97t has to be replaced by 0j, lc i = 
e/ (— uoC + Gi + G r ), where we denote the capacitance of the system by C. The modification to the 
current is analogous to Eq. (I200> . 

Finally, as far as the experiments are concerned, a In T behavior in the linear resistivity of short metal- 
lic bridges, together with an I- ^-characteristic which agrees well with the universal function dl90> has 
been observed in II 131 . The bridges have been about L — 80 nm long, 80 nm wide and 10 nm thick; 
the temperature region where the logT in the linear resistivity has been seen is between T w 100 mK 
and T w 2 K. In 111 131 the Coulomb correction to the tunneling conductance has been suggested as the 
explanation; the considerations presented here and in| 1 12 1 show that, also in the intermediate regime with 
both diffusive and interface resistivity the predicted /-^-characteristic does not change and thus agrees 
with the experimentally observed one. In [ 1 1 3 1 the Thouless energy, which sets the scale for the lowest 
diffusive mode in an open system and therefore for the low temperature saturation of the conductance, has 
been estimated to be of the order of several Kelvin, whereas the lnT is observed down to 100 mK. In 
the case with resistive interfaces, on the other hand, the energy of the lowest diffusive mode is reduced, 
7o ~ HD / L 2 (R wlrc / R) <C hD/L 2 . Although some caution is necessary in applying these considerations 
to the experiment 2 , we conclude that in the measurement of 1 1 13 1 the diffusive resistance may be consider- 
ably smaller than the interface resistance. A further hint for the importance of interfaces is found from the 
prefactor A: For the open system we found A as 0.25, and A — 2G\G r j (G r + G r ) 2 in the tunnel limit. 
The experimental values 1 1 13 1 are between A w 0.43 . . . 0.7, i.e. closer to the tunnel limit than to the open 
system. 



5.4 Diffusive point contacts 

In the previous subsection we observed that the electron escape rate 1 / r csc = 70 determines the temper- 
ature below which the Coulomb interaction contribution to the conductivity saturates. The reason for this 
saturation is the assumption that the system is connected, on both sides, to "ideal leads" where quantum 
fluctuations are completely suppressed. In this subsection we will allow for voltage fluctuations in the 
leads, i.e. fluctuations of the field </>j nt (x, t), which lead to a temperature dependent conductance also in 
the low temperature region ksT < h/r csc . When the system under consideration consists of a tunnel 
junction, the fluctuations in the leads are dominated by the conventional Coulomb interaction effects on 
the tunnel conductance as we discussed in subsection l5.ll 

However, an interesting question is whether a tunnel junction is strictly necessary see Coulomb blockade 
physics lll4llll5llll6lfTT7llll8l . As far as transport properties are concerned the issue was addressed 
by Golubev and Zaikin 1 1 17 1, and Levy Yeyati et al. II 181 . who studied within a scattering approach the 
transport through a coherent scatterer coupled to an electromagnetic environment. For a single conducting 
channel, the correction to the conductance due to the interaction with environmental modes was found to 



The interface resistance required to obtain the low temperature saturation of the conductance at ~ 100 mK is about one order 
of magnitude larger than the actually measured resistance; on the other hand there is still room to refine the theory. 
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be (T = 0) 

2e 2 ,„ , f , RcZ _ 2e 2 mn „ 

— = — r-(l-r / dw , G=— r, (201) 

G n J eV ijj h 

where r is the channel transmission and the impedance Z characterizes the environment. For a poorly 
transmitting channel (r <C 1) this agrees with standard Coulomb blockade theory in the perturbative 
limit, compare Eqs. ( !157l i-( fT6ll . Apparently in a highly transmitting channel the relative change in the 
conductance is suppressed by the factor (1 — r) with respect to a tunnel junction. In the multichannel case 
the reduction factor (1 — r) is replaced by II 171 

(3= En T n(l-T n ) ^ (2Q2) 



The parameter (3 is known from the theory of shot noise 1 123 1; it is equal to one for tunnel junctions, and 
(3 = 1/3 for diffusive conductors. 

We will now derive the correction to the current through a diffusive conductor due to the interactions in 
the leads within the quasiclassical theory. When one neglects the Coulomb interaction inside the wire, the 
change of the Green's functions due to the interaction in the leads fulfills the equations 

{-yuj + Ddl) 5g z {r,;* ,u) = (203) 

(-iw - Ddl)5g K {r r ^,uj) = 0. (204) 



The corresponding solutions decay or increase exponentially on the scale = y/ ' D/ui. When the struc- 
ture under consideration is longer than L w , the fluctuations in the leads influence 8g Z ' K only in a narrow 
region near the interfaces, and will only weakly affect the transport. In the opposite limit, when the struc- 
ture is shorter than L u or when considering frequencies below the Thouless energy D/L 2 , the approximate 
solutions of the equations are the linear interpolation of g z and g K between the values on the left and right 
interfaces, i.e. 

6g z (i T ,x,t) » ^—^Sg z (r 1 ;0,t) + jSg z ( m L,t) (205) 

V%x,i) » £—l S g K ( m 0,t) + ^5g z ( m L,t). (206) 

One observes that the Coulomb correction to the current through the wire can be expressed in terms of the 
correction to the Green's functions at the interfaces between the wire and the leads. In the following, for 
simplicity, we consider interactions just in one of the two leads. Under this condition we determine the 
correction to the current by inserting Eqs. ( 1205 \ and d206i into the general expression for diffusive systems 
obtained in Eq. dl05> . and integrate between and L, with the result 

5I JSEl = ~G f de^M[F(e,L)-F(e,0)], (207) 
6e J N Q 

where G is the conductance of the system in the absence of fluctuations in the leads. The result obtained 
within the quasiclassical method is consistent with the scattering approach mentioned above: Similar to 
what is found for a tunnel junction, the Coulomb interaction correction to the conductance is related to 
the tunneling density of states at the interfaces also for a fully coherent diffusive metal. For a diffusive 
conductor the variation of the conductance due to the interaction is reduced by a factor three compared to 
tunnel junctions. 

The density of states at the interface is sensitive to the geometry of the leads. When the leads have a 
three-, two-, or one-dimensional character, the density of states in the leads is proportional to 

yft d = 3 

SN(e) oc { ln(e) d = 2 (208) 
lA/e d=l 
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at zero temperature. Notice, however, that the density of states on the surface of a three-, two-, or one- 
dimensional conductor is not identical to the density of states in the bulk. For example, for a semi-infinite 
three-dimensional sample the density of states as a function of energy and of the distance from a metal- 
insulator interface is 1 103 1 



6N 3d (e,z) = ^ o 1 +_ !_ ( 2 P— e-^cosz + lnter)) (209) 
y/2 2n 2 D 3 / 2 8n 2 Dz \ J XQ x '} 



with x = z^2e/D, x\ = z^/2/Dt. The energy dependent part of the above equation varies thus by a 
factor two between the surface and the bulk, 



5N? d 



surface 



^=2^(6) = ^^372, (210) 

and the typical length scale for this variation is L e — y/D/e. We expect that 5./V s 3 ^ rfacc (e) as it is given here 
is also applicable to an interface between a three-dimensional lead and a diffusive point contact, provided 
the diameter of the point contact is small compared to L t . 

In order to discuss the one-dimensional case we first go back to the general expression for the density 
of states correction, 

SN(e,x) = N Im / ^ / dxi.F(e - w;xi)p (x, cj;xi)</> lnd (xi, u;x). (211) 

In metals the field 0i n d usually propagates much faster than the charge density po\ in particular, in the limit 
of good screening, the charge neutrality condition corresponds to a field 0i n( j which propagates instantly 
over the full system under consideration. As a consequence the field 0; n d(xi, u>; x) is in many situations 
a smooth function of the position, compared to the charge density po(x, u>; xi) which is concentrated in 
a region |x — xi I - y/D/u). In one dimension, for example, the induced potential varies linearly in the 
position with the typical scale of the system size, see for example Eq. d!94t . When neglecting the spatial 
variation of both 0i n d and the distribution function, then the variable xi in the equation above can easily 
be integrated to yield the density of states as 

<W(e,x) « N lm /"j£ F ( c _ w;x )_l_ ^ nd ( x ,w;x). (212) 
J 2ir —no 

When the leads are in thermal equilibrium the first and second derivative of the current with respect to the 
voltage are (T = 0) 

§ = G+ i G «>, (213) 

ov 3 A^o 

and 

d 2 I _ 2G ReZ( e y,0) 

dV 2 3 V h/e 2 ' 1 ; 

In the last equation we replaced the induced electrical potential by the local impedance, </>; n d(x, u>; x) = 
eZ(cu, x). Remember that 0i n d(x, t; x) is the voltage at (x, t), due to a charge that has been placed at x and 
t = 0. At the interface between a point contact and leads the voltage decays due to the charge flowing into 
the leads or into the contact, and therefore Z(oj, x = 0) is the total impedance of the leads and the point 
contact in parallel. For example, when we describe the point contact by a capacitance C and a resistance 
R, and the leads as a i?C-transmission line with a resistance Rq and a capacitance Co per unit length, the 
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impedance Z is given by 1 94 1 



Z 



1 



•-1 

lead 



+ z, 



'-1 

point contact 



(215) 



z 



r-1 

point contact 



-\wC + R~ 



(216) 



R L if L < y/l/uC Ro 



(217) 



v /i?o/ - iwC ifL> y/l/u,C Ro 



For long leads, where Z(lo) is proportional to 1/y/uj, the density of states correction acquires its usual 
one-dimensional form, SN(e) oc 1/y/e. 

The situation where Zi ea d(a;) is frequency independent has been analyzed in detail by Golubev and 
Zaikin in 111 171 . with the result 



with R~{ = Rr 1 + {RqL)- 1 , and 7 « 0.577. 

Golubev and Zaikin suggested that the In T behavior in the conductance, and the zero bias anomaly 
observed by Weber et al. 11 131 in short metallic bridges, might be due to the effects described above; in 
this case one concludes that the electromagnetic environment of the bridges behaves as an Ohmic resistor, 
with the environmental resistance at least comparable to or larger than the resistance of the nanobridge 
(R ~ 10 il); for related experiments see also 1 1 19 1. 

In concluding this section we summarize the basic mechanisms by which the Coulomb interaction 
affects transport through small structures. At low temperatures - in the absence of inelastic scattering - 
the interaction effects may be grouped in two classes. First, the internal electric potential created in a 
diffusive conductor leads to an additional source of random scattering that interferes with scattering from 
impurities. This is a quantum effect and depends on the details of charge diffusion. We have discussed in 
subsections !5 ,2l and l5 . 3l how charge diffusion is modified in small structures, and how this in turn modifies 
transport. Second, adding a charge to a conductor costs energy and therefore the transport is blocked. This 
is avoided in macroscopic samples via screening, but the latter can be inefficient in small structures, where 
the necessary flow of the background charge is slow due to the geometric constraints. We have discussed 
this type of physics in subsections 15 . 1 1 and 15 .41 where the screening properties were taken into account in 
terms of an effective environmental impedance Z. 

Of course, in real systems, the two groups of interaction effects may be difficult to separate. However, 
in a large number of experimentally relevant cases, this separation is indeed possible due to the different 
time scales for screening and for charge diffusion. In ordinary metals, for instance, screening is much 
faster than charge diffusion. 

6 Summary 

In this article we discussed transport properties of mesoscopic samples with emphasis on electron-electron 
interaction effects. The theory was formulated in terms of the quasiclassical Green's functions. These 
Green's functions are averaged over impurities from the beginning, and therefore the complex diagram- 
matic analysis is avoided. We demonstrated that many known effects can be obtained within this formalism: 
the interaction corrections to the diffusive transport, density of states, and also non-perturbative Coulomb 
blockade physics. 

A major part of this article has been devoted to the regime of diffusive transport, where we worked 
out a theory for electron-electron interaction effects far from thermal equilibrium. There is a long list of 
possible extensions: for example it would be interesting to study time dependent phenomena. Response 
to time dependent fields with frequencies in the microwave regime has been studied in some detail for 
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weak localization, both theoretically Ml 3I IT20 30 1 and experimentally 1121111221 . In these experiments the 
influence of a microwave field on the DC-conductivity of thin films was investigated. In the experimental 
studies, besides weak localization, the interaction effect also is present, but there are only few theoretical 
studies of the interaction correction in the presence of high frequency fields. Altshuler et al. | 13 1 calculated 
the high-frequency (linear) conductivity. A first attempt to calculate the DC-conductivity in the presence 
of microwaves was made by Raimondi et al. 1271 . There the microwave field was included in the theory 
via the diffusion propagators. The distribution function, on the other hand, was considered stationary. 
Going beyond this analysis, the distribution function could be determined by explicitly solving the kinetic 
equation in the presence of the microwave. The thus determined distribution function has then to be used 
in the expression for the interaction correction to the current density, Eq. M 1 31 . 

In Section[5]we discussed transport properties of simple devices such as films and wires with interfaces. 
This immediately suggests that our formalism may be applied to hybrid systems, where also interfaces 
between different materials are present, e.g. metal-superconductor, metal-ferromagnet, etc. Furthermore 
we restricted ourselves to the limit of diffusive electron motion; in small and clean samples, however, the 
mean free path can be longer than the systems size, which requires a theory for ballistic electron motion. 
A theory for the Coulomb correction to transport in ballistic systems has recently been put forward by 
Zala et al. 1531 in order to describe transport at intermediate temperatures, where the thermal time tt is 
comparable to the (elastic) scattering time r. At present, however, the theory has not yet been applied to 
finite size systems. 

We conclude the list of possible extensions of the work presented here by noting that we have concen- 
trated on the ensemble averaged current density, leaving aside the rapidly developing field of current noise 



In relation to the problem of electron dephasing, we have shown that dephasing affects also the interac- 
tion correction to the conductivity. However there are still many open questions: What is the microscopic 
nature of the low temperature saturation of the dephasing time? What is the microscopic origin of the 
dephasing rate which varies as \[T in gold films 1731 . and as T 1 / 3 in E^Si^CuOg [77 1? Why is the 
dephasing rate, as determined from the weak localization magneto-resistance, different from the rate de- 
termined from conductance fluctuations measurements 1701 1711 ? 

To summarize, transport in mesoscopic- and nano-structures is an active field of research, with several 
open questions ranging from fundamental problems of quantum theory to possible technological applica- 
tion in the future. 

Acknowledgements We acknowledge discussions with C. Castellani, T. Luck, and U. Eckern. This work was sup- 
ported by the German Italian exchange program Vigoni (CRUI and DAAD) and miur-cofin2002022 534. 

A Field theoretic approach: the non-linear sigma model 

The field theoretic formulation of the interacting, disordered electron system was pioneered by Finkelstein 
in the 1980's |46 1. Recently, this formulation has been extended to the non-equilibrium case by means of 
the Keldysh technique by Kamenev and Andreev [51 1, and Chamon, Ludwig and Nayak [52 1 in the case of 
normal-conducting metals, and by Feigel'man, Larkin and Skvortsov 1 125 1 in the case of superconductors. 
Given the already extensive literature available on the subject, we believe that, rather than repeating the 
derivation of the non-linear sigma model, it is perhaps more useful, to show how the formalism applied in 
the present paper - namely quasiclassical Green's function technique in the presence of fluctuating internal 
fields - is found when starting from the non-linear sigma model. 

In the absence of the electron interaction, the action of the non-linear sigma model is given by 



H23lll()6lfT24l . 
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where the trace is a summation or integration over all degrees of freedom, and the time derivative is here 
understood to be dtQ = dtQttx \ti=t- The field Q = Q^ t , (x) must satisfy the constraint Q 2 = 1, i.e. 



E 



dt Qf tl (x)Q* , (x) =S(t- t')8 ik . (220) 



An external scalar field cxt and the Hubbard-Stratonovich field $ due to the Coulomb interaction are 
described by the following terms in the action: 

iSi = -i7riVoTr$ Q 7 Q Q + Tr$ Q l/-VS /3 $ /3 (221) 
\S 2 = -i7riV o Tre0 cxt Q + 47VoTre0 cxt $2, (222) 

with 7 1 = (T° (the 2x2 unit matrix) and j 2 = <r x , and V is the "statically screened" interaction, i.e. a 
constant. 

In order to make contact with the quasiclassical approach, we express the Green functions in terms of 
the Q-fields. First we observe that 1511 



G tt -(x,x') = 



i -i 



Go" 1 + ttQ + $ «7 a 
It 



>«,», (223) 



where the brackets (...) q $ indicate that one has to average over the fields Q and $. By using the condition 
that the fields Q and <f> are only slowly varying in space and time, one finds a relation between the £- 
integrated Green function and the Q-matrix. In particular, upon taking the s-wave component one finds 



(224) 



At this point we make the crucial approximations: At first we restrict the field Q to a ^-dependent 
saddle point configuration. Varying the action with respect to Q, under the condition Q 2 = 1, we find 



1 



irN Q Dd 2 Q - irN d t 8(t - <i)£(x - xi) - i7r7V (e<; 



$«7 a )+AQ + QA = 0. 



(225) 



where A = X\ J tl (x) is a Lagrange multiplier. Making use of the normalization condition and after elimi- 
nating A the saddle point equation reads 



d t Q - Dd x (Qd x Q) - i [e^ ext + $ q7 q , Q] = 0, 



(226) 



which agrees with the kinetic equation for the quasiclassical Green's function in the presence of an internal 
field. Hence we identify Q with g. In the next step we restrict $ to a saddle point plus quadratic fluctuations. 
Variation of $) with respect to $ gives the saddle point for the internal field, 



$f F (x,i) = VttNo 



Q t 1 t 2 (x)--e0 cxt (x,i) 







(227) 
(228) 



In order to obtain this relation we used that Q 21 = 0, and Q\} + Q 22 = at the saddle point. The field 
fluctuations are determined from the second derivative of the action with respect $, 



d 2 S 

d$2 



d SS\6Q SSS 2 Q 



d SS 



d$ SQ J <5$ SQ 5<I> 2 d$ 5$ ' 



(229) 
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Since 5S/SQ = for all field configurations $ the first two terms in the equation above are equal to zero. 
The final result is 



-i(^x,i)$,(x',t')} = -V£ yn (x,t;x',0, 
where the dynamically screened interaction V dyn is 



[y^]- 1 = y- 1 5(x-x')^-o( J I 



+ 



x fl (x,t;x',0 
X A (x,t;x',t') X *(x,*;x',0 



with 



X K (x,t;x',t') = -7dV 



X A (x',t';x,i) 



£0i(x',t') 



(230) 

(231) 

(232) 
(233) 



x *'(x,,;x', t <) = . 

<50 2 (x',i') d02(x',t') 

B Fermi liquid parameters and the spin triplet channels 

In the entire article we took into account the electron spin only by factors of two appearing at various 
places. When considering the electron spin explicitly, the Green's function acquires spin indices, and the 
Zeeman term has to be taken into account in the equation of motion: 



d t g - Dd x (gd x g) - i[e<fi + [i B cr ■ B, g] = 0. 



(234) 



Both the scalar field <fi and the magnetic field B are a sum of internal and external contributions, <fi = 
4> cyit 1 + int , B = B ext 1 + B int . By separating the spin singlet and triplet components of the Green's 
function as 

g t 



g5 ss , +g- er ss >, g = \-'^g ss , g = ^ ^ er ss ,g s , s 



2 

s s 

the charge density and magnetization density are expressed as 

p(x,i) = 2eN 



- 5 12 (x,t)-e0 1 (x, t) 



m(x, t) = 2/j, B N 



,12 



(x, t) - n B B 1 (-x,t) 



(235) 

(236) 
(237) 



where (fix, B x are the sum of the external field and the diagonal component of the corresponding internal 
field. As in the spinless case we assume that the internal field has Gaussian fluctuations around some mean 
value. The mean field equations read 



6 MF 



MF 



= An 



AO, 



IT 1 

2 e 



7T 1 
2 (IB 



g t Y(x)-B c 



(238) 
(239) 



On average, the quantum components of the fields are zero, (fi\ lF = B| IF = 0. In the equations above Aq 
and Aq are the Landau parameters, which have the following relations to the standard parameters Fq' s , 
and to the interaction constants in the spin singlet and triplet channels, 7 s t : 
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Before discussing the fluctuations we demonstrate how to obtain the charge susceptibility within the mean 
field scheme. The variation of the charge density in the presence of an applied external field i 



is 



2eN n 



,12 



(x)-e(l 



(x,i)-^-^(x) 



(242) 
(243) 



here we made use of the self-consistency condition for the internal field, Eqs. j238t and ( I239l l. The variation 
of the Keldysh component of the Green's function due to the applied field is determined by the equation of 
motion. At first we note that the commutator of the field with the Keldysh Green's function at equal times 
is 



2i 



{[d>,g 12 ]) tt = --d t m- 
By using this relation it is found that 



(-kj + Dq 2 )5g 12 (q, u ) = -(-ho) [e0 oxt (q, w) + e< F (q, w)] , 

and then 

[-Ml - A$) + Dq 2 } <fo 12 (q,w) = - A^e^^u;). 

Finally, we find the charge susceptibility in the presence of Fermi liquid renormalizations, 

Dq 2 



X (q,o;) = 2e 2 iVo(l-^)- 



- A s ) + Dq 2 



(244) 



(245) 



(246) 



(247) 



In particular the static susceptibility is modified by the factor (1 — Aq) in the presence of the interac- 
tions. Similar results are obtained for the spin susceptibility. The susceptibility, obtained here using the 
quasiclassical formalism, agrees with the direct diagrammatic approach, see e.g. 1471 . 

The field fluctuations are also fully determined within our formalism, as we discussed in the previous 
section. Here we give the expression in the absence of an external magnetic field, and for a spatially 
homogeneous system. The retarded interaction in the spin singlet channel is 



»2 = 



-ho + Dq 2 



N -i(l - A s )oj + Dq 2 
and in the spin triplet channel, 

A^ -ho + Dq 2 



No -i(l - Ag)w + Dq 



(248) 



(249) 
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